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Some  approaches  are  developed  to  analyze  and  interpret  wide  angle  bottom  loss 
data  in  terms  of  ocean  sediment  acoustic  properties.  After  a  brief  review  of 
bottom  loss  data  collection  and  the  state-of-the-art  knowledge  of  relevant 
ocean  bottom  properties,  a  filter  description  of  the  acoustic  signal  is 
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time  versus  grazing  angle  plots.  It  is  shown  that  these  effects  can  be 
identified  and  that,  in  general,  they  severely  distort  the  ocean  bottom 
impulse  response.  -  •  ,~J 
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20.  (Cont'd) 

The  autocorrelation  and  power  cepstrum  are  examined  for  their  usefulness  in 
extracting  geoacoustic  information  from  the  multiple  bottom  bounce  data. 

The  autocorrelation  is  initially  preferred  due  to  several  reasons  having 
to  do  with  the  fact  that  it  has  the  superposition  property  with  respect  to 
convolution: 

(1)  the  shape  of  events  are  more  interpretable  in  terms  of 
the  phase  of  the  time  components, 

(2)  the  amplitudes  of  events  are  more  predictable,  and 

(3)  the  effects  of  coherent  noise  are  more  recognizable. 

On  the  other  hand,  the  power  cepstrum  is  as  good  and  better  than  phe 
autocorrelation  for  much  synthetic  plane  wave  data,  since  it  gives  better 
resolution  of  delays  between  components  (it  produces  a  spike  at  the  delays 
between  components  rather  than  a  pulse  with  the  crosspower  spectrum  of  the 
components)  and  since  it  is  better  able  to  take  advantage  of  infinite 
dynamic  range.  The  advantage  of  infinite  dynamic  range  is  lost  and  the  power 
cepstrum  is  degraded  with  respect  to  the  autocorrelation  by  adding  white 
uncorrelated  noise  to  synthetics.  For  real  data,  preference  of  one  technique 
over  the  other  will  be  dependent  on  source-receiver  geometry,  geoacoustic 
structure,  and  signal-to-noise  ratio. 

An  impulse  response  estimate  was  obtained  from  the  data  by  using  a  least 
squares  method  to  deconvolve  an  (n-1)  bottom  bounce  waveform  from  an  (n) 
bottom  bounce  waveform. n  This  method,  which  is  directly  analogous  to  the 
debubbling  technique  of  Wood  et  al.  (1978),  requires  that  the  two  waveforms 
correspond  to  equal  source  and  receiver  depths  and  equal  bottom  grazing 
angles.  Although  good  results  have  been  obtained,  significant  noise  contami¬ 
nation  of  the  deconvolved  signal  results  for  source  depth  difference  between 
3  and  7%.  This  noise  contamination  can  be  reduced  by  low  pass  filtering  at 
the  cost  of  resolution. 

Application  of  the  autocorrelation,  power  cepstrum,  and  deconvolution 
techniques  to  real  data  obtained  in  a  turbidite  province  indicated  that  bottom 
loss  in  the  region  was  largely  controlled  by  a  single  reflector  and  refracted 
energy.  Although  ambiguity  exists  as  to  whether  the  reflector  is  a  bottom  or 
subbott<im  reflector,  either  interpretation  implies  a  significant  difference 
of  geoacoustic  structure  from  that  previously  predicted  for  this  region. 
Assuming  the  two  components  of  the  impulse  response  are  the  specularly 
reflected  and^refracted  components,  the  gradient  for  the  first  180  m  decreases 
from  2.0  sec  to  1.3  sec-*  with  an  average  of  about  1.5  sec-*  in  the  first 
180  m.  The  attenuation  coefficient  was  about  0.02  dB/m-kHz  in  this  same  depth 
interval.  The  autocorrelations  and  power  cepstrums  indicated  ificant 
reflectivity  in  the  upper  50  m  of  the  subbottom. 
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CHAPTER  1 


BACKGROUND 


1 . 1  Acoustic  Propagation  Experiments 

The  purpose  of  this  work  is  to  develop  approaches  to  analyze  and 
interpret  wide  angle  multiple  bottom  bounce  acoustic  signals  in  terms  of 
ocean  sediment  acoustic  properties.  This  falls  within  the  realm  of  the 
science  of  underwater  sound  which  has  developed  as  a  quantitative  subject 
in  the  years  since  World  War  II  (Urick,  1975).  Recently  in  underwater 
acoustics  interest  has  been  focused  on  long  range  propagation  in  deep 
water.  Often  these  acoustic  environments  are  bottom  limited  so  that  a 
significant  fraction  of  the  acoustic  energy  interacts  with  the  ocean  bot¬ 
tom  (Hampton,  1978).  It  is  desired  to  characterize  the  physics  of  this 
interaction  with  the  ocean  bottom  and,  in  particular,  to  be  able  to  predict 
the  amount  of  loss  expected  at  the  ocean  bottom  boundary  as  a  function  of 
angle  and  frequency.  This  involves  a  quantitative  knowledge  of  geoacous¬ 
tic  structure  as  a  function  of  depth. 

Experiments  that  have  been  designed  to  study  acoustic 
propagation  in  bottom  limited  environments  are  known  as  bottom  loss  exper¬ 
iments.  These  experiments  are  generally  long  range  (1-200  km),  low  fre¬ 
quency  (0-600  Hz)  measurements  with  a  geometry  like  that  shown  in  Fig.  1. 
Vertical  arrays  are  deployed  with  hydrophone  clusters  at  approximately  1400 
and  1900  m  below  the  surface  and  30  m  above  the  bottom,  with  one  hydrophone 
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located  directly  on  the  bottom.  The  arrays  are  anchored  to  the  bottom  and 
suspended  by  a  float. 

Explosive  SUS  (shallow  underwater  source)  charges  are  then 
commonly  deployed  by  a  ship  moving  away  from  the  receiver  site  and  deto¬ 
nated  at  depths  of  approximately  18,  91,  and  24*1  m.  Figure  1  shows  that 
for  the  depths  of  most  ocean  basins  (3.5-5  km)  most  of  the  acoustic  signals 
measured  in  the  experiment  have  interacted  with  the  bottom  several  times 
before  reaching  the  receiver  (Mitchell  et  al.,  1978). 

In  this  work  it  is  desired  to  make  quantitative  inferences  of  the 
acoustic  properties  of  ocean  bottom  sediments  from  these  acoustic  measure¬ 
ments. 

1 .2  Bottom  Loss  Analysis  Procedure 

From  a  theoretical  standpoint,  the  interaction  of  acoustic 
signals  with  the  ocean  bottom  can  be  described  by  a  single  function,  the 
reflection  coefficient,  R(5,w),  a  function  of  grazing  angle  and  frequency. 
This  function  is  commonly  calculated  using  a  normal  mode  model  (Hawker  and 
Foreman,  1978)  or  a  ray  tracing  model  (Mitchell  and  Lemmon,  1979)  by 
considering  a  single  frequency  infinite  plane  wave  interacting  with  the 
ocean  bottom. 

The  ocean  bottom  is  modeled  as  a  series  of  flat,  horizontal 
(generally  fluid)  layers  with  the  variables  of  layer  thickness,  density, 
sound  speed,  sound  attenuation,  and  sound  speed  gradient.  If  the  sediments 
have  significant  rigidity  then  shear  wave  parameters  versus  depth  must  be 
included  in  the  model  (Vidmar  and  Foreman,  1978;  Vidmar,  1980).  Bottom 
loss  (BL)  is  defined  as  negative  twenty  times  the  common  logarithm  of  the 
plane  wave  reflection  coefficient: 
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-20  log  R(0,u>) 


where  R(0,w)  is  the  reflection  coefficient. 

The  measurement  of  bottom  loss  within  the  constraints  of  this 
definition  is  virtually  impossible.  In  practice  bottom  loss  is  obtained 
from  propagation  loss  experiments  in  which  the  acoustic  wave  has  reflected 
from  the  ocean  bottom  one  or  more  times.  Figure  2  is  a  flowchart  which 
shows  the  process  of  deriving  bottom  loss  from  multiple  bottom  bounce 
acoustic  data  (Mitchell  et  al.,  1978). 

First,  the  recorded  signals  are  digitized  and  stored  on  magnetic 
tape.  Then  the  approximate  arrival  timing  structure  of  waterborne  paths  is 
calculated  from  ancillary  data  using  nominal  source  and  receiver  depths 
and  the  water  sound  speed  structure.  Third,  the  arrivals  are  plotted  and 
the  predicted  arrival  structure  verified.  Finally,  the  spectra  of  the 
reference  structure  are  calculated,  corrected  for  geometrical  spreading 
loss,  and  compared  with  the  spectra  of  data  to  obtain  bottom  loss.  The 
source  spectra  levels  for  the  data  discussed  in  this  report  and  processed 
at  ARL:UT  were  assumed  to  be  those  reported  by  Gaspin  and  Shuler  (1971). 
The  ocean  bottom  is  assumed  to  be  flat  and  horizontal  so  the  grazing  angle 
at  the  ocean  bottom  is  unique  and  can  be  calculated  from  the  source- 
receiver  geometry  and  the  water  sound  speed  structure. 

Figure  3  is  a  plot  of  an  intermediate  step  of  the  ARL:UT 
processing  of  bottom  loss  with  the  results  of  other  steps  summarized.  The 
digitized  multiple  bottom  bounce  signal  is  plotted  above  the  time  axis  with 
two  tick  marks  equaling  1  sec.  The  top  printed  line  documents  the  number 
of  bottom  bounces  for  each  blip,  which  was  determined  on  the  basis  of 
travel  time  information.  It  is  seen  that  arrivals  with  less  than  two 
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bottom  bounces  are  essentially  canceled  due  to  simultaneous  arrival  of 
ghosts  at  these  times.  The  grazing  angle  for  each  blip  is  calculated  from 
the  sound  speed  structure  in  the  water  column  and  the  source-receiver  geom¬ 
etry.  The  reference  structure  is  calculated  from  assumed  source  levels  and 
the  calculated  arrival  structure.  The  grazing  angle  for  each  blip  and  the 
reference  loss  (at  50  Hz)  for  each  blip  is  tabulated  in  the  figure.  The 
reference  loss  includes  the  effects  of  geometrical  spreading  of  the  acous¬ 
tic  energy  and  absorption  by  the  seawater.  The  measured  loss  is  calculated 
by  windowing  each  blip  and  finding  the  difference  between  the  reference 
level  and  the  measured  level.  Bottom  loss  is  simply  the  difference  between 
the  reference  and  measured  losses.  The  measured  loss  and  bottom  loss  (both 
at  50  Hz)  are  tabulated  for  the  signal  of  Fig.  3  below  the  time  axis. 
Generally,  bottom  loss  is  calculated  at  frequencies  of  25,  50,  100,  200, 
and  400  Hz  with  a  1/3  octave  band  average.  This  gives  bottom  loss  as  a 
function  of  angle  and  frequency. 

To  obtain  source  levels,  some  bottom  loss  experiments  use 
measured  levels  of  direct  arrivals  that  have  not  interacted  with  the  ocean 
bottom  rather  than  assuming  levels  from  previous  work  (Santaniello  et  al., 
1979).  This  involves  adding  the  predicted  spreading  loss  of  the  direct 
arrival  (WPLp)  to  the  level  of  the  direct  arrival  (RLp) 

SL  =  RLd+WPLd  .  (1.2.1) 
Bottom  loss  (BL)  can  then  be  calculated  by  subtracting  the  level  of  the 
bottom  arrival  (RLg)  and  its  predicted  spreading  loss  (WPLg)  from  the 
source  level. 

BL  =  SL  -  RLg  -  WPLg  ,  (1.2.2) 
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or 

BL  =  (RLp-RLg)  +  (WPLd-WPLb)  (1.2.3) 
Expression  (1.2.2)  is  the  expression  used  at  ARL:UT  to  derive  bottom  loss 
using  the  source  levels  predicted  by  Gaspin  and  Shuler.  Alternatively, 
Santaniello  et  al.,  with  an  experimental  geometry  designed  to  allow  for 
indirect  measurement  of  the  source  levels,  use  expression  (1.2.3)  to  cal¬ 
culate  bottom  loss. 

For  multiple  bottom  bounce  r".  *ls  bottom  loss  must  be 
normalized  to  the  loss  per  bottom  bo># ..he  bottom  loss  calculated 
from  (1.2.2)  or  (1.2.3)  for  an  (n)  bot'1  n  'n.-.-.nce  signal  must  be  divided  by 

n. 

1 . 3  Reflected  and  Refracted  Energy  in,  the  Ocean  Bottom 

In  order  to  predict  the  amount  of  energy  loss  expected  at  the 
ocean  bottom  in  different  ocean  basins  a  physical  model  is  needed  to 
predict  the  main  energy  components  returned  to  the  water  column.  The 
simplest  model  is  one  that  considers  the  ocean  bottom  as  a  perfect 
absorber.  In  this  case,  all  of  the  energy  which  interacts  with  the  ocean 
bottom  is  absorbed  by  it.  Alternatively,  the  ocean  bottom  can  be  modeled 
as  a  reflector  with  a  simple  Rayleigh  reflection  coefficient;  all  energy 
not  reflected  at  the  water-sediment  interface  is  assumed  to  be  absorbed 
into  the  ocean  bottom.  These  models  yield  white  transfer  functions  for  the 
ocean  bottom. 

A  more  complex  model  is  to  consider  a  refracted  component  which 
returns  energy  to  the  water  column  due  to  a  positive  sound  speed  gradient. 
This  results  in  a  two-component  model  of  the  ocean  bottom  transfer  function 
with  one  component  reflected  at  the  water-sediment  interface  and  another 
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component  transmitted  at  the  water-sediroent  interface  and  refracted  back 
to  the  water  column  through  the  sediments  (see  Figs.  4  and  5).  The  exis¬ 
tence  and  the  importance  of  this  refracted  arrival  has  been  known  for  many 
years . 

Apparently,  Hill  (1952)  was  the  first  to  identify  an  arrival 
refracted  through  the  ocean  subbottom.  His  interpretation  was  later  sup¬ 
ported  by  a  more  controlled  experiment  conducted  by  Officer  (1955)  and  by 
interpretations  of  bottom  loss  experiments  of  Morris  (1970)  and  Hanna 
(1973).  Christensen  et  al.  (1975)  conducted  an  experiment  in  which  they 
were  able  to  characterize  the  refracted  arrivals  in  terms  of  travel  times, 
amplitudes,  and  frequency  content. 

Since  that  time  controlled  experiments  employing  deconvolution 
processing  (Dicus,  1976;  Herstein  et  al.,  1979;  Santaniello  et  al.,  1979) 
have  revealed  that  the  refraction  arrival  is  a  significant  event  espe¬ 
cially  for  low  frequencies.  Santaniello  et  al.  attribute  the  calculation 
of  negative  bottom  loss  to  the  breakdown  of  the  plane  wave  assumption  and 
the  occurrence  of  refracted  energy.  The  amplitude  of  the  refracted  arrival 
i3  frequency  dependent  due  primarily  to  the  frequency  dependent  attenua¬ 
tion  of  the  sediments.  The  specular  reflection,  however,  should  contain 
more  high  frequencies  and  may  be  nearly  white.  The  overall  transfer 
function  of  this  model  will  be  colored. 

The  two-component  model  described  above  assumes  that  the 
impedance  in  the  subbottom  is  a  smooth  and  continuous  function  of  depth; 
however,  discontinuous  or  sharp  impedance  contrasts  with  depth  yield 
reflectors  in  the  subbottom  besides  the  specular  reflection.  Many  recent 
studies  indicate  that  there  are  significant  subbottom  reflectors  in 
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addition  to  the  specularly  reflected  and  refracted  component.  For 
instance,  in  the  work  of  Herstein  et  al.  in  the  Hatteras  Abyssal  Plain 
(1979)  and  Santaniello  et  al.  (1979)  at  least  one  significant  subbottom 
reflector  was  observed.  The  data  presented  by  Dicus  (1976;  1981)  also 
showed  significant  subbottom  reflectivity.  Data  obtained  in  the  Tufts 
Abyssal  Plain  presented  by  Chapman  (1980)  indicated  that  the  bottom  loss 
was  controlled  by  two  major  reflectors;  refracted  energy  did  not  appear  on 
the  deconvolved  traces.  Data  of  Santaniello  et  al.  (1979)  also  indicate 
that  subbottom  reflected  energy  dominates  the  bottom  return  rather  than 
refracted  energy.  It  appears  then  that  the  model  of  the  ocean  bottom  as  a 
specular  reflection  with  a  refracted  component  is  not  accurate  for  some 
ocean  environments,  and  knowledge  of  sharp  and/or  discontinuous  impedance 
contrasts  with  depth  may  be  required  to  predict  bottom  reflection  loss  in  a 
particular  ocean  basin. 

A  primary  motivation  of  the  work  presented  here  is  the  desire  to 
examine  subbottom  reflectivity  for  different  ocean  basins  using  existing 
bottom  loss  data.  Currently,  the  two  component  model  of  a  specularly 
reflected  and  refracted  arrival  is  often  U3ed  to  match  measured  bottom  loss 
to  theoretical  curves  using  geoacoustic  models  without  sharp  subbottom 
impedance  contrasts.  In  what  follows,  two  correlation  techniques  and  a 
deconvolution  technique  are  applied  to  multiple  bottom  bounce  signals 
obtained  in  a  turbidite  province.  These  techniques  reveal  two  major  compo¬ 
nents;  however,  it  is  ambiguous  whether  the  reflected  component  is  the 
specular  reflection  or  a  subbottom  reflection.  Either  interpretation 
implies  that  the  geoacoustic  structure  in  this  turbidite  province  is 
significantly  different  from  that  proposed  a  priori  (Hamilton,  1979). 
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1.4  Geoacoustic  Structure 


There  are  four  main  parameters  in  the  ocean  bottom  which  affect 
bottom  loss:  the  depth  and  nature  of  the  basement  (roughness,  impedance, 
etc.)  sound  speed  versus  depth,  density  versus  depth,  and  attenuation  ver¬ 
sus  depth.  The  most  important  of  these  parameters  is  the  depth  and  nature 
of  the  basement.  If,  for  instance,  the  basement  is  shallow  (<75  m)  and 
rough,  as  in  much  of  the  open  ocean  where  the  basement  is  basaltic  and 
undulatory  (Fig.  6),  then  the  acoustic  response  of  the  basement  will 
largely  control  the  level  and  shape  of  the  bottom  loss  curves.  It  would  be 
expected  that  the  frequency  dependence  of  the  acoustic  response  of  the 
basement  will  be  highly  dependent  on  the  spatial  frequency  distribution  of 
the  basement  roughness  and  that  the  dependence  of  bottom  loss  on  grazing 
angle  will  be  dependent  largely  on  the  depth  of  the  basement  below  the 
ocean  bottom.  In  the  case  where  the  sediment  thickness  is  large  (>75  m) 
and  the  basement  is  smooth  and  hard,  then  the  frequency  response  of  the 
basement  could  be  considered  to  be  white,  and  its  effect  on  bottom  losr 
curves  would  be  seen  mainly  at  low  frequencies  and  high  grazing  angles. 

The  most  important  parameter  after  the  depth  and  nature  of  the 
basement  is  the  sound  speed  versus  depth.  The  sound  speed  versus  depth 
controls  the  "paths"  that  acoustic  energy  travels  between  source  and 
receiver  and,  therefore,  the  densities  and  attenuations  it  encounters. 
Thus,  the  dependence  of  bottom  loss  on  the  grazing  angle  of  the  incident 
acoustic  wave  is  very  dependent  on  the  sound  speed  structure  in  the 
subbottom. 
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At  low  grazing  angles  (0-10°)  bottom  loss  is  highly  dependent  on 

the  sound  speed  ratio  at  the  water-sediment  interface  (c  /c  ).  This  is  due 

w  s 

partly  to  the  intromission  effect  where  (for  a  sound  speed  ratio  less  than 
one)  the  specular  reflection  amplitude  goes  to  zero  at  some  small  grazing 
angle.  More  importantly  there  is  an  interference  effect  between  the  specu¬ 
larly  reflected  and  refracted  arrivals.  This  interference  effect  gives  a 
large  peak  in  the  bottom  loss  curve  for  sound  speed  ratios  less  than  one  at 
small  grazing  angles.  Bottom  loss  goes  to  zero  as  the  grazing  angle  goes 
to  zero  for  sound  speed  ratios  greater  than  one. 

The  gradient  of  sound  speed  versus  depth  mainly  controls  the 
level  of  the  bottom  loss  curves  for  grazing  angles  not  associated  with 
basement  interaction.  This  is  due  to  the  fact  that  the  gradient  of  sound 
speed  controls  the  path  lengths  of  the  transmitted  paths  and  the  attenua¬ 
tion  and  density  contrasts  that  the  acoustic  energy  encounters.  Disconti¬ 
nuities  in  velocity  versus  depth  (or  sharp  changes)  will  result  in  reflec¬ 
tions  that  also  will  significantly  affect  the  level  and  shape  of  the  bottom 
loss  curves  for  angles  not  associated  with  bottom  interaction. 

Attenuation  is  the  least  well  known  of  the  geoacoustic 
parameters  versus  depth,  and  this  is  due  partly  to  the  uncertainty  of  the 
relative  importance  of  different  mechanisms  which  cause  acoustic  energy  to 
be  lost  in  the  ocean  bottom.  It  is,  however,  a  very  significant  parameter 
in  predicting  bottom  loss  at  all  grazing  angles. 

The  importance  of  the  density  profile  on  the  shape  and  level  of 
the  bottom  loss  curves  is  important  at  discontinuities  or  sharp  changes 
with  depth  which  result  in  reflections  such  as  at  the  water-sediment 
interface. 
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Edwin  Hamilton,  who  has  been  greatly  involved  in  geoacoustic 
modeling,  has  recently  given  a  lengthy  review  article  of  the  methodology 
and  state-of-the-art  knowledge  of  geoacoustic  modeling  of  the  ocean  floor 
(Hamilton,  1980).  Figure  7  is  a  set  of  velocity  profiles  presented  by 
Hamilton  for  different  sediment  types  from  reflection  and  refraction  pro¬ 
filing.  This  figure  shows  the  approximate  variability  of  velocity  versus 
depth  in  the  ocean  floor.  It  is  seen  that  the  smallest  increase  in 
velocity  with  depth  occurs  with  siliceous  sediments,  the  largest  increase 
occurs  with  calcareous  sediments,  and  an  intermediate  increase  occurs  with 
terrigenous  turbidites. 

These  profiles  are  obtained  primarily  from  seismic  reflection 
and  refraction  profiles.  The  actual  velocity  profile  for  an  ocean  basin  is 
determined  generally  by  linear  regression  of  interval  velocities  obtained 
from  seismic  reflection  data,  which  is  forced  through  a  sound  speed  value 
at  the  water-sediment  interface  obtained  from  physical  property  measure¬ 
ments  of  geologic  cores. 

Figure  8  shows  velocity-density  data  indicating  the  relationship 
between  velocity  and  density  for  different  sediment  types.  In  geoacoustic 
modeling,  density  versus  depth  is  often  determined  by  obtaining  a  velocity 
profile  and  converting  the  velocities  to  densities  using  velocity-density 
relationships  in  conjunction  with  known  or  assumed  sediment  types. 

Figures  9-11  show  published  predictions  of  attenuation  versus 
depth.  Figure  9  is  the  composite  of  laboratory  measurements,  physical 
property  measurements  of  cores,  and  other  field  data  collected  from  acous¬ 
tic  measurements.  Note  that  attenuation  varies  over  a  considerable  range 
of  0.015  dB/m-kHz  to  0.35  dB/m-kHz.  Jacobson  et  al.  report  a  similar  range 
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FIGURE  7 

INSTANTANEOUS  VELOCITY  versus  DEPTH  CURVES 
FOR  DIFFERENT  SEDIMENT  TYPES 


[Reproduced  from  Hamilton,  (1980)1 
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FIGURE  10 

ATTENUATION  versus  DEPTH  IN  THE  BAY  OF  BENGAL 

[Reproduced  from  Jacobson  et  al.,  (1981)] 
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FIGURE  11 

ATTENUATION  AS  A  FUNCTION  OF  DEPTH  IN  THE 
SEA  FLOOR  CALCULATED  FROM  BOTTOM  LOSS  DATA 

(Reproduced  from  Mitchell  and  Focke,  (1980)1 
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of  values  from  a  seismic  refraction  experiment  in  the  Bay  of  Bengal,  in 
good  agreement  with  Hamilton's  data  (Fig.  10).  On  the  other’  hand  Fig.  11 
shows  attenuation  values  obtained  from  bottom  loss  data  b'  Mitchell  and 
Focke  (1980)  that  are  considerably  lower  than  those  obtained  by  Jacobson  or 
Hamilton.  At  present  there  is  no  explanation  for  the  difference  in  the  two 
profiles.  It  could  be  due  to  actual  differences  in  intrinsic  attenuations 
of  the  sediments  in  the  two  regions  or  to  use  of  erroneous  assumptions  in 
determining  the  attenuation  values.  The  majority  of  the  current  litera¬ 
ture  seems  to  point  to  the  lower  values  of  attenuation  at  least  for  shallow 
(<300  m)  depths  (Bottom  Interacting  Acoustics.  198O). 

Figures  12  and  13  illustrate  the  probable  need  for  refinement  of 
preliminary  models  of  velocity  and  density  versus  depth.  This  is  because 
preliminary  geoacoustic  models  often  do  not  take  into  account  sharp  imped¬ 
ance  contrasts,  and  even  more  rarely  account  for  lateral  variability  at  a 
site.  Figures  12  and  13  show  data  obtained  from  physical  property 
measurements  of  the  Deep  Sea  Drilling  Project  (DSDP)  for  cores  obtained  in 
the  Gulf  of  Mexico.  Hole  90  was  drilled  in  the  western  Gulf  and  Hole  91 
was  drilled  about  200  km  away  in  the  northern  Gulf.  Figure  12  shows 
density  measurements  obtained  by  the  DSDP  and  Fig.  13  shows  velocity  data 
calculated  from  the  density  measurements  using  the  velocity-density  rela¬ 
tionships  of  Hamilton  and  the  sediment  types  reported  from  the  cores. 

Assuming  that  the  variance  of  the  errors  of  measurement  is 
constant  between  sites,  it  is  seen  that  the  velocities  are  higher  in  the 
northern  Gulf  than  the  western  Gulf  and  that  there  are  more  sharp  impedance 
contrasts.  The  substantial  agreement  between  the  velocity  profile  derived 
from  seismic  measurements  in  the  western  Gulf  with  the  DSDP  measurements 
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gives  credence  to  both  and  implies  that  errors  in  the  DSDP  velocity  values 
obtained  from  cores  are  small.  Thus,  these  data  illustrate  significant 
lateral  variability  that  occurs  over  a  fairly  small  areal  extent.  These 
data,  particularly  Fig.  13(a),  also  imply  the  feasibility  of  comparing 
acoustic  experiments  and  coring  data  and  suggest  the  use  of  DSDP  data  for 
preliminary  geoacoustic  modeling. 

Finally,  Figs.  14  and  15  show  the  effects  of  different 
geoacoustic  environments  on  averaged  bottom  loss  curves.  Figure  14  shows 
bottom  loss  for  an  area  of  thin  sediment  overlying  a  rough  basaltic  base¬ 
ment.  It  is  seen  that  there  is  relatively  high  loss  at  low  angles  and  that 
the  loss  does  not  increase  with  frequency  for  all  grazing  angles. 

Figure  15  shows  two  bottom  loss  plots  for  thick  sediment 
regions.  It  is  seen  that  the  loss  increases  with  frequency  for  all  grazing 
angles  and  is  lower  at  the  low  grazing  angles  than  the  thin  sediment 
region.  It  is  seen  that  there  is  higher  loss  for  the  calcareous  ooze 
region  than  for  the  region  of  terrigenous  turbidites.  The  reason  for  this 
difference  is  ambiguous.  It  could  be  due  to  a  higher  velocity  gradient  in 
the  turbidite  region,  a  higher  reflectivity,  or  lower  attenuations. 

The  work  that  follows  is  an  attempt  to  develop  processing 
techniques  that  could  resolve  the  ambiguity  of  the  bottom  loss  curves  by 
identifying  events  of  the  time  signal  and  thereby  constraining  the  model  of 
the  geoacoustic  structure. 

1 .5  Forward  Modeling 

In  this  study  it  was  often  necessary  to  model  acoustic  signals  as 
a  function  of  time  and  grazing  angle.  These  synthetic  waveforms  are  gen¬ 
erated  from  a  computer  program  called  IMPULSE  (and  its  feeder  program 
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FIGURE  14 

BOTTOM  LOSS  IN  AN  AREA  OF  THIN  SEDIMENT 
OVERLYING  LOUGH  BASEMENT 
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FIGURE  15 

BOTTOM  LOSS  CURVES  FOR  TWO  AREAS  OF  THICK  SEDIMENT 


BOTREF),  which  resulted  from  the  work  of  Mitchell  and  Lemmon.  (See  "A  ray 
theory  model  of  acoustic  bottom  interaction,"  J.  Acoust.  Soc.  Am.  66 
(1979).)  IMPULSE  calculates  the  synthetic  waveforms  by  calculating  a  syn¬ 
thetic  source  waveform  by  the  method  of  Gaspin  and  Schuler  (197“'),  adding 
in  ghosts  at  the  respective  ghost  delays,  and  convolving  this  "total  source 
waveform"  with  the  plane  wave  ocean  bottom  impulse  response  calculated  by 
the  ray  tracing  model  described  below. 

1.5.1  Impulse  Response  Calculation 

The  impulse  response  of  the  synthetic  waveforms  discussed  in 
this  thesis  were  calculated  by  tracing  rays  between  layers  of  the  subbottom 
and  calculating  reflection  and  transmission  coefficients  at  each  interface 
in  the  subbottom  with  the  usual  fluid-fluid  and  fluid-solid  boundary  con¬ 
ditions  (Mitchell  and  Lemmon,  1979).  The  ocean  bottom  is  modeled  as  an 
isovelocity  fluid  half-space  (water)  above  an  arbitrary  number  of 
absorbing  fluid  layers  with  variable  sound  velocity,  density,  and  attenua¬ 
tion  which  overlie  a  lower,  solid  half-space  (basement)  (see  Fig.  17(a)). 

A  simple  example  of  the  method  is  shown  in  Fig.  16.  In  this  case 
an  isovelocity  half-space  overlies  an  absorbing  fluid  half-space  with  a 
positive  velocity  gradient..  It  is  seen  that  there  is  a  reflected  compo¬ 
nent,  a  refracted  component,  and  multiple  refracted  components  which  have 
reflected  from  beneath  the  interface  one  or  more  times.  At  any  point  along 
the  interface  there  is  an  incident  ray,  a  reflected  ray,  an  upcoming  ray, 
and  a  downgoing  ray.  The  field  can  be  represented  by  the  displacement 
potentials,  <p j,  <j>R,  4>D,  and  d>y.  Realizing  that  the  downgoing  wave  at  point 
A  is  the  same  as  the  downgoing  wave  at  0  to  within  a  phase  factor,  taking 
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(a)  GEOACOUSTIC  MODELING  FOR  INPUT  TO  BOTREF 


(b)  RAY  DIAGRAM  ILLUSTRATING  METHOD  OF 
IMPULSE  RESPONSE  CALCULATION  OF  BOTREF 


FIGURE  17 

IMPULSE  RESPONSE  CALCULATION  OF  BOTREF 

[Reproduced  from  Mitchell  and  Lemmon,  (1979)1 
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into  account  absorption,  and  the  t/2  phase  shift  associated  with  a  plane 
wave  passing  through  a  turning  point  (Silbiger,  1968;  Tolstoy,  1 968 ; 
Burash,  1968),  the  reflected  field  can  be  found  to  be: 

P2(1  +  a)k1sinei  -  P1  ( 1-a)k2sin02 
R  “  p£(  1  +  ‘»)k1sin^1  +  ( 1-a)k2sin02  ’ 


where 

P 1  2  z  density  in  upper  (lower)  interface 
k.  _  =  wavenumber  in  upper  (lower  interface) 

'  t 

ff.  _  =  incident  (transmitted)  angle,  and 

*  f  ^ 

is  the  factor  which  takes  into  account  the  phase  shifts  and  attenuation 
of  the  refracted  arrival  relative  to  the  reflected  arrival. 

«  =  exp  (-  /  B(z)dz  -  i»(rA0-rB0>  -  i^) 

Mitchell  and  Lemmon  then  generalizes  this  result  to  a  series  of  flat 
horizontal  layers  by  tracing  rays  and  solving  the  fluid-fluid  and  fluid- 
solid  boundary  conditions  for  a  plane  wave  locally  at  each  interface  (see 
Fig.  17(b)).  The  results  of  this  procedure  are  the  recursive  relations. 
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where 


C  ..sine 

-  (I,  ,,  /a  -ay  R  ,  .  r--~-r  y  - 

(n,n+1)  n  n  n  n,n+1)  C  smj  , 

n  n-1 


In  n+^  =  incident  amplitude, 

Rn  n+1  =  re^lecte<^  amplitude, 

Cn  ?n  =  velocity  and  density  at  the  top  of  the  nth  interface,  and 

c  .p  ,  =  velocity  and  density  at  the  bottom  of  the  (n-1) 
n- I  n- i 

interface, 

and  where  a  ,  y  are  terms  which  take  into  account  attenuation  and  time 
n  n 

delays  between  impulse  response  components. 


Thus,  the  plane  wave  reflection  coefficient  of  the  stack  of  layers  is 


Figures  18  and  19  show  the  results  of  the  Mitchell  and  Lemmon  work. 
Figure  18(a)  shows  the  Fourier  transform  of  the  reflection  coefficient  at 
a  particular  angle  of  incidence.  This  is  the  transfer  function  of  the 
ocean  bottom.  Figure  18(b)  shows  the  bottom  loss  prediction  from  BOTREF  as 
a  function  of  angle  and  frequency.  Finally,  Fig.  19  shows  how  this  ray 
tracing  model  compares  to  normal  mode  theory.  Figure  19(a)  shows  that  the 
ray  tracing  model  is  good  for  frequencies  greater  than  about  10  Hz. 
Figure  19(b)  shows  the  breakdown  of  the  model  when  there  is  a  critical 
angle.  The  advantage  of  the  ray  tracing  model  over  the  normal  mode  model 
is  that  the  impulse  response  can  be  calculated  much  more  rapidly  on  the 


computer 
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BOTTOM  LOSS  -  dB 


(a)  REFLECTION  COEFFICIENT  FOR  A  TURBIDITE 
MODEL  AT  20°  GRAZING  ANGLE 


(b)  BOTTOM  LOSS  AT  50  Hz  FOR  A  TURBIDITE 
MODEL  WITH  A  5°  CRITICAL  ANGLE 


FIGURE  19 

COMPARISON  OF  WAVE  THEORY 
AND  RAY  THEORY  PREDICTIONS 
OF  BOTTOM  RESPONSE 

[Reproduced  from  Mitchell  and  Lemmon,  (1979)] 
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In  what  follows,  the  synthetic  ocean  bottom  transfer  function 
will  be  modeled  by  the  ray  tracing  method.  Figure  20  shows  the  impulse 
response  output  from  BOTREF  using  the  geoacoustic  model  of  Table  I  over 
300  msec  with  the  specular  reflection  aligned  at  25  msec.  In  the  absence 
of  sharp  impedance  constrasts  with  depth,  these  specular  reflected  and 
refracted  components  are  the  only  important  components  of  the  ocean  bottom 
response,  and  knowledge  of  the  time  difference  between  these  components  is 
sufficient  to  obtain  a  velocity  profile  with  depth. 

1.5.2  Synthetic  Source  Waveform 

The  source  waveform  that  is  used  in  synthetic  modeling  of 
acoustic  signals  in  this  report  comes  from  the  work  of  Gaspin  and  Shuler 
(1971),  who  model  the  source  waveforms  of  explosive  charges  by  fitting 
exponentials  to  the  amplitudes  of  positive  and  negative  peaks  of  the  source 
signature.  The  Fourier  transform  of  the  constructed  waveform  provides  an 
estimate  of  the  source  spectrum  levels  for  bottom  loss  processing,  and  the 
constructed  waveform  itself  is  used  in  the  synthetic  modeling  of  acoustic 
signals.  The  shock  pulse  is  modeled  as  an  exponential  pulse  of  the  form 


where  PQ  is  the  amplitude  of  the  shock  pulse  and  J  is  a  decay  constant 
evaluated  from  empirical  relationships.  The  bubl'  e  pulses  are  modeled  as 
double  exponential  decays. 

The  parameters  involved  in  the  formulation  of  the  synthetic 
waveform  are  shown  in  Fig.  21.  It  is  seen  that  if  the  source  waveform  is 
considered  to  have  four  significant  bubble  pulses,  there  are  five 
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amplitude  measurements,  four  bubble  pulse  delays,  and  three  periods 
associated  with  the  width  of  the  bubble  pulses  that  must  be  determined 
empirically.  The  delays  between  bubble  pulses  were  found  by  Gaspin  and 
Shuler  to  be  proportional  to  the  cube  root  of  the  charge  weight  and 
inversely  proportional  to  the  sixth  root  of  the  depth  of  detonation. 

Although  this  provides  a  practical  and  convenient  method  for  the 
estimation  of  the  source  spectrum,  it  should  be  noted  that  although  an 
actual  explosion  shock  wave  decays  exponentially  the  decay  rate  decreases 
markedly  after  a  time  of  about  2.  The  shock  wave  has  a  significant  amount 
of  energy  in  the  negative  phase  with  the  time  of  the  first  zero  crossing 
strongly  dependent  on  the  depth  of  burst.  The  bubble  pulses  are  also  not 
well  modeled  by  exponentials. 

Thus,  although  these  waveforms  are  useful  for  synthetic  studies 
in  which  primarily  the  bubble  pulse  timing  and  amplitude  must  be  modeled 
adequately,  they  are  not  accurate  enough  models  to  use  in  source  deconvolu¬ 
tion  of  real  data. 

1.5.3  Synthesized  Waveform 

To  derive  the  synthesized  waveform,  the  effects  of  ghosts  must 

» 

be  added  to  the  source  waveform.  This  "total  source  waveform"  must  then  be 
convolved  with  the  ocean  bottom  impulse  response  to  obtain  the  synthesized 
waveform. 

Ghost  time  delays  can  be  calculated  by  ray  tracing  using  the 
experiment  geometry  and  the  water  sound  speed  structure.  Assuming  that  the 
transfer  function  of  the  air-water  interface  is  a  negative  one,  the  "total 
source  waveform"  is  calculated  by  adding  the  source  waveforms  at  zero  delay 
and  each  ghost  delay  with  the  appropriate  sign,  positive  or  negative. 
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With  the  additional  assumption  that  the  grazing  angle  at  the 
ocean  bottom  is  the  same  for  the  first  arrival  and  its  ghost  arrivals,  the 
synthesized  waveform  is  calculated  by  convolving  the  "total  source  wave¬ 
form"  with  the  impulse  response.  Since  the  bottom  is  assumed  to  be  flat 
and  laterally  homogeneous,  the  impulse  response  for  an  n  bottom  bounce 
signal  will  be  the  Fourier  transform  of  (R(0,o>))n  (-1)n"t'1  where  R(0,w)  is 
the  ocean  bottom  transfer  function  for  a  one  bottom  bounce  signal. 

The  assumption  that  the  impulse  response  for  the  first  arrival 
and  its  ghost  arivals  is  the  same  becomes  more  valid  as  the  number  of 
bottom  interactions  increases. 


CHAPTER  2 


*» 


TEM  FILTER  EFFECTS  ON  ACOUSTIC  DATA 


2.1  System  Filter  Effects 

If  the  impulse  response  consists  essentially  of  a  specularly 
reflected  and  a  refracted  arrival,  knowledge  of  the  travel  time  difference 
between  these  components  would  allow  for  the  calculation  of  velocity  as  a 
function  of  depth  in  the  subbottom.  Measuring  this  travel  time  difference 
from  bottom  loss  data  is  difficult,  for  common  source  receiver  geometries 
of  the  experiment,  because  source  bubble  pulses  and  the  source  and  receiver 
ghosts  considerably  obscure  the  time  delay  of  interest.  In  this  section 
the  arrivals  of  the  source  and  receiver  ghosts  and  source  bubble  pulses 
will  be  shown  on  time  versus  grazing  angle  plots  for  multiple  bottom  bounce 
acoustic  signals. 

Figure  22  shows  the  frequency  representation  of  an  acoustic 
signal  recorded  in  a  bottom  limited  ocean  and  Fig.  23  shows  a  schematic 
that  illustrates  the  various  filtering  effects  which  act  on  the  signal.  It 
is  seen  that  the  signal  consists  of  three  filter  components:  the  source 
waveform,  the  source  and  receiver  ghosts,  and  the  multiple  bottom  bounce 
impulse  response. 

First,  consider  the  filtering  effect  of  source  bubble  pulses  on 
time  versus  grazing  angle  plots.  Since  the  bubble  pulse  delays  are  depend¬ 
ent  on  source  characteristics  only,  the  effect  of  the  source  bubble  pulses 
will  be  independent  of  grazing  angle  or  range.  Thus,  the  interfering 
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FREQUENCY  DOMAIN  REPRESENTATION 
OF  ACOUSTIC  SIGNAL 


ARL:UT 
AS-81325 
Gl  GA 
3  4  81 


OCEAN  BOTTOM 


FIGURE  23 

SYSTEM  FILTERING  EFFECTS 


ARL 


effects  due  to  bubble  pulses  of  aligned  arrivals  will  appear  as  straight 
lines  on  time  versus  grazing  angle  plots  at  the  bubble  pulse  delays.  In 
fact  every  event  of  the  impulse  response  will  have  a  trail  of  bubble 
pulses. 

For  the  work  done  here  the  nominal  source  depths  are  18.3,  91.4, 
and  244.1  m  (or  60,  300,  and  800  ft).  Figure  24  shows  the  quasi- 
theoretical  plots  of  these  waveforms  for  the  first  300  msec  from  the  work 
of  Gaspin  and  Shuler.  Figure  20  shows  that  we  are  attempting  to  measure 
time  delays  in  the  impulse  response  which  vary  from  0  to  about  125  msec  for 
grazing  angles  between  0  and  30°.  Recognizing  the  convolutional  nature 
of  filtering  operations  in  the  time  domain,  it  is  seen  that  only  the 
18  m  source  will  not  considerably  obscure  the  time  delay  of  interest. 

Second,  we  consider  the  filtering  effects  due  to  ghosts.  From 
Fig.  23  it  is  seen  that,  for  each  acoustic  pulse  generated  at  the 
source,  four  arivals  come  in  at  the  receiver:  the  arrival  that  is  not 
reflected  from  the  air-water  interface,  the  source  ghost,  the  receiver 
ghost,  and  the  source  plus  receiver  ghost. 

The  delays  of  these  ghosts  can  be  approximated  by  considering  a 
plane  wave  interacting,  with  a  free  surface  (see  Fig.  25).  This  time 
delay  is 

.  _  2d  sing 
tg  =  c 

where  d  is  the  depth  of  the  source  or  receiver,  c  is  the  water  sound 
speed,  and  9  is  the  grazing  angle. 

This  delay  is  illustrated  on  a  time  versus  grazing  angle  plot 
in  Fig.  26  for  depths  between  12.5  and  400  m.  Zero  delay  is  aligned  at 
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FIGURE  26 


25  msec  to  be  consistent  with  plots  of  real  and  synthetic  acoustic  signals 
to  be  shown  in  later  sections. 

Figure  27  shows  the  acoustic  response  of  the  ghosting  effects. 
The  source  ghost  response  is 

6(t)  -  5(t-ts) 

and  similarly  the  receiver  ghost  response  is 

6(t)  -  *(t-tR) 

The  convolution  of  these  filter  components  yields  the  total  ghost  response 

6(t)  -  i(t-ts)  -  «(t-tR)  +  a(t-tg-tR) 

which  consists  of  the  source  ghost,  the  receiver  ghost,  and  the  source  plus 
receiver  ghost.  As  expected  the  source  ghost  and  the  receiver  ghost,  which 
reflect  from  the  air-water  interface  only  once,  are  negative  while  the 
other  two  arrivals  are  positive.  The  total  ghost  response  shows  that 
significant  interference  will  occur  at  the  source  ghost  time  delay,  the 
receiver  ghost  time  delay,  and  their  sum.  Since  the  time  delay  is  propor¬ 
tional  to  depth,  for  the  plane  wave  approximation,  significant  interfering 
energy  will  arrive  along  the  lines  drawn  in  Fig.  26  at  the  source  and 
receiver  depths  and  the  sum  of  their  depths. 

The  final  filter  needed  to  describe  multiple  bottom  bounce 
acoustic  signals  is  the  multiple  bottom  bounce  impulse  response.  The  (n) 
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bottom  bounce  impulse  response  is  equal  to  the  convolution  of  the  (n-1) 
bottom  b.^ijce  impulse  response  and  the  negative  of  the  one  bottom  bounce 
impulse  response  (assuming  the  bottom  is  flat  and  laterally  homogeneous). 
If  the  ocean  bottom  impulse  response  consists  of  only  a  specular  reflection 
and  a  refracted  component  which  has  gone  through  a  caustic,  then  it  can  be 
represented  as  a  pulse  followed  by  a  pulse  with  a  90°  phase  shift.  This  is 
the  case  illustrated  in  Fig.  28  (BB  =  bottom  bounce).  It  is  seen  that  in 
this  case  there  is  a  90°  phase  shift  between  each  of  the  components  of  the 
multiple  bottom  bounce  impulse  response.  The  fact  that  the  leading  arrival 
is  alternately  positive  and  negative  is  due  to  the  reflection  at  the  air- 
water  interface  between  bounces.  Thus,  the  phase  of  the  impulse  response 
components  of  multiple  bottom  bounce  signals  can  be  predicted. 

Amplitudes  can  also  be  predicted.  Consider  the  same  case  of  an 
impulse  response  that  consists  essentially  of  a  specularly  reflected  and 
refracted  component.  If  we  consider  the  delay  between  the  components  to  be 
unity  the  signal  can  be  represented  in  the  Z  domain  by 

.  .  ,  .  iTr/2Z  , 

+  a2U’)e 

where  a1  is  the  amplitude  of  the  specular  reflector,  a2  is  the  amplitude 
of  the  refracted  arrival,  and  Z  is  the  unit  delay  between  them. 
Dividing  by  a1  and  recognizing  that  convolution  in  the  time  domain  is 
equivalent  to  multiplication  in  the  Z  domain  we  have: 
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where  a=^a2(uj)/a1(j)j  e1  "/2.  Thus,  it  is  seen  for  the  case  where  a2/a.|  =  1 
that  the  impulse  response  will  be  symmetric  and  have  its  greatest 
contribution  in  the  center.  If  a,,/a^>1  then  the  impulse  response  will 
have  its  larger  contributions  in  its  tail. 

So  it  is  seen  that  the  amplitude  and  phase  of  the  components  of 
the  impulse  response  can  be  predicted  for  multiple  bottom  bounce  signals 
if  the  one  bottom  bounce  impulse  response  is  known.  The  fact  that  the 
phase  of  the  impulse  response  components  can  be  predicted  gives 
encouragement  that  these  components  can  be  identified  in  the  data.  This 
would  allow  for  the  inversion  of  the  time  delay  between  them  to  velocity 
and  attenuation  profiles  as  a  function  of  depth  in  the  subbottom. 

To  get  an  idea  of  how  the  multiple  bottom  bounce  impulse 
response  components  are  delayed  as  a  function  of  grazing  angle,  consider 
a  model  of  the  ocean  bottom  as  an  infinite  half-space  with  a  simple 
linear  gradient  of  velocity.  In  this  case  the  travel  time  of  the 
refracted  path  in  the  sediment  is 

t  .  I,  I  c°3gR  \ 

TR  "  "  g  ln  y  1  +  sin0Rj  ’ 

where  g(sec_1)  is  the  gradient  and  0R  is  the  transmitted  angle  at  the 
interface  (Ewing  and  Leet,  1932).  In  the  plane  wave  case,  assuming  the 
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water  column  is  locally  isovelocity  at  the  water-sediment  interface,  the 
travel  time  difference  is 


At 


qo3*R  \ 
1  +  sinflp  ) 
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-  sine 
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(Mitchell  and  Lemmon,  1980.)  In  reality,  the  spreading  of  acoustic  energy 
from  the  source  is  more  nearly  spherical  than  planar.  The  main  effect  of 
this  is  that  the  reflected  and  refracted  components  have  different  angles 
of  incidence  at  the  water-sediment  interface.  The  plane  wave  assumption 
will  be  used  since  it  pr  ides  a  good  estimate  of  the  delay  between  the 
reflected  and  refracted  components  and  is  very  easily  calculated  from  the 
preceding  formula.  Appendix  I  shows  a  comparison  between  the  plane  wave 
and  a  very  simple  spherical  wave  model  of  At  versus  grazing  angle. 

The  plane  wave  assumption  is  made  in  all  of  the  work  of  the  next 

section. 

2.2  T- 0  Plots  of  Multiple  Bottom  Bounce  Acoustic  Signals 

In  this  section  the  predicted  travel  time  curves  of  the 
components  of  multiple  bottom  bounce  acoustic  signals  will  be  synthesized 
and  compared  to  measured  signals  with  the  first  recognizable  phase  of  the 
measured  signals  aligned  at  25  msec.  This  will  illustrate  the  effects  of 
the  source  waveform  and  ghosts  on  travel  time  versus  grazing  angle  plots  of 
acoustic  signals,  and  will  show  that  these  effects  can  often  be  identified 
and  that,  in  general,  they  severely  distort  the  ocean  bottom  impulse 
response.  (The  aligned  arrival  is  not  usually  believed  to  be  the  first 
arrival  of  the  acoustic  signal  since  the  first  arrival  is  usually  buried  in 
the  noise;  rather  it  is  the  first  phase  that  can  easily  be  recognized  on 


all  traces.  Aligning  an  easily  recognized  phase  is  time  consuming  ('- 2  h). 
Aligning  less  easily  recognized  phases  would  be  much  more  time  consuming 
and  the  alignment  of  a  particular  phase  would  be  less  certain.) 
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Figure  29  shows  the  predicted  travel  time  curves  for  the  arrivals 
of  a  five  bottom  bounce  impulse  response  assuming  a  plane  wave  model  with  a 
water-sediment  sound  speed  ratio  of  0.994  and  a  gradient  of  1.1  sec-1. 
This  is  a  good  first  order  prediction  of  the  travel  time  curves  of  the 
bottom  response  components  based  on  geoacoustic  models  proposed  in  the 
literature.  The  important  components  of  the  impulse  response  are  assumed 
to  be  a  single  specular  reflection  and  a  single  refracted  component  per 
bottom  bounce,  as  shown  in  Fig.  28.  Clearly,  a  more  complicated  model  with 
subbottom  reflectors  will  considerably  complicate  the  travel  time  plot  and 
make  unambiguous  identification  of  arrivals  more  difficult.  Internally 
reflected  events  (or  deep  refracted  events)  are  generally  ignored  since 
they  are  considered  to  be  attenuated  to  the  noise  level. 

Figure  30  shows  the  effects  of  the  source  bubble  pulses  of  the 
previous  figure  with  a  244  m  source.  For  all  pass  plots  (0-600  Hz)  the 
shock  pulse  would  be  expected  to  be  largest  with  an  exponential  decay  of 
the  amplitude  of  trailing  bubble  pulses.  For  a  low  pass  plot  the  first 
bubble  pulse  may  be,  and  often  is,  observed  to  be  larger  than  the  shock 
pulse.  In  general,  the  second  bubble  pulse  is  very  small  and  the  third 
bubble  pulse  can  generally  be  ignored. 

Figure  31  shows  the  ghosting  effects  of  the  multiple  bottom 
bounce  impulse  response  on  a  travel  time  versus  grazing  angle  plot.  Each 
component  of  the  impulse  response  is  delayed  by  the  ghost  delays  associated 
with  the  source  depth,  the  receiver  depth,  and  the  sum  of  the  source  and 
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FIGURE  30 

BUBBLE  PULSE  DELAYS  OF  MULTIPLE 
BOTTOM  BOUNCE  IMPULSE  RESPONSE 

(244  m  SOURCE) 


receiver  depths.  The  depths  of  the  source  and  receiver  were  chosen  in 
accordance  with  the  nominal  depths  of  the  data  that  follow.  Figure  32 
shows  the  combination  of  all  filtering  effects:  the  source  bubble  pulses, 
ghosting  effects,  and  the  multiple  bottom  bounce  impulse  response. 

Figure  33  shows  data  with  nominal  source  and  receiver  depths  of 
244  and  498  m,  respectively.  A  particular  phase  is  aligned  at  25  msec. 
Figure  34  shows  that  the  aligned  phase  is  probably  not  the  first  component 
of  the  acoustic  signal.  A  precursor  of  the  aligned  arrival  is  observed 
with  variable  advance.  The  significance  of  this  is  that  not  all  of  the 
signal  components  will  be  delayed  with  respect  to  the  aligned  arrival  at 
25  msec. 

Also  in  Fig.  34  a  phase  is  identified  that  is  delayed  with  respect 
to  the  peak  aligned  at  75  msec.  This  time  delay  seems  to  be  a  mirror  image 
of  the  time  delay  of  the  precursor  identified  in  Fig.  33.  This  gives 
credence  to  the  existence  of  impulse  response  components  with  this  time 
difference.  Thus,  delays  between  impulse  response  components  apparently 
can  be  identified  on  aligned,  otherwise  unprocessed,  acoustic  signals. 
However,  this  time  difference  is  more  easily  identified  and  interpreted 
using  techniques  which  follow. 

Returning  to  Fig.  33  it  is  seen  that  there  are  three  trends  of  the 
acoustic  energy  on  the  time  versus  grazing  angle  plot.  The  first  trend  is 
located  on  a  line  rotated  at  about  10°  with  respect  to  the  vertical;  the 
second  trend  is  located  on  a  line  rotated  at  about  35°  with  respect  to  the 
vertical,  and  the  third  trend  is  located  on  a  line  rotated  about  45°  with 
respect  to  the  vertical.  The  first  trend  is  due  to  refracted  energy  of  the 
impulse  response,  while  the  other  two  trends  are  due  to  ghost  energy. 
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FIGURE  32 

TRAVEL  TIME  versus  GRAZING  ANGLE  PLOT 
OF  TIME  DELAYS  BETWEEN  MAIN  ARRIVALS 
OF  A  5  BB  ACOUSTIC  SIGNAL  WITH  A  244  m 
SOURCE  AND  498  m  RECEIVER 
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GRAZING  ANGLE -deg 


REFRACTED 


FIGURE  34 

SIX  BOTTOM  BOUNCE  ACOUSTIC  SIGNALS  WITH  244  m  SOURCE  DEPTH 
AND  498  m  NOMINAL  RECEIVER  DEPTH  ALIGNED  AT  75  msec  TO 
REVEAL  AT  LEAST  ONE  ARRIVAL  BEFORE  THE  ALIGNED  ARRIVAL 
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Figure  35  shows  the  6  BB  data  of  Figs.  33  and  34  aligned  at 

25  msec  with  travel  time  curves  associated  with  the  five  bottom  ice 

impulse  response  (c  /c  =0.994,  g  =  1.1  sec”1),  the  source  waveform,  and 
s  w 

the  source  and  receiver  ghosts.  The  5  BB  impulse  response  was  used  rather 
than  the  6  BB  response  since  at  least  one  arrival  has  been  tentatively 
identified  before  the  aligned  arrival.  Thus,  at  least  one  component  of  the 
impulse  response  will  not  be  delayed  with  respect  to  the  aligned  arrival. 
The  nominal  source  and  receiver  depths  were  used  to  predict  bubble  pulse 
delays  and  source  and  receiver  ghost  delays. 

From  Fig.  35  it  is  is  seen  that  there  are  agreements  and 
disagreements  between  the  predicted  travel  times  and  the  data.  In  partic¬ 
ular,  the  vertical  lines  representing  the  bubble  pulse  delays  of  the 
aligned  arrival  agree  with  the  bubble  pulse  periodicity  seen  in  the  data. 
Also,  in  agreement,  the  travel  time  curves  for  the  source  ghost  overlie  a 
significant  trend  of  energy  in  the  data.  On  the  other  hand,  in  disagree¬ 
ment,  there  is  a  large  amount  of  acoustic  energy  predicted  in  the  lower 
right  hand  corner  of  the  plot  where  the  data  show  none.  Also,  in  disagree¬ 
ment,  there  is  a  large  arrival  at  about  90  msec  at  30  °  grazing  angle  that 
is  not  predicted  by  the  travel  time  curves.  The  first  disagreement  is  due 
to  the  fact  that  the  nominal  receiver  depth  is  in  error.  The  second 
disagreement  is  probably  due  to  the  fact  that  the  a  priori  choice  of 
gradient  is  in  error. 

Figure  36  shows  data  recorded  at  the  same  receiver,  but  due  to  an 
18  m  source.  The  advantage  of  the  18  m  data  is  in  the  identification  of 
filter  components  such  as  the  receiver  ghost  and  the  bubble  pulses.  Thus, 
with  the  18  m  source  data  it  will  be  easier  to  show  that  the  nominal 
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FIGURE  36 

FIVE  BOTTOM  BOUNCE  ACOUSTIC  SIGNALS 
WITH  18  m  SOURCE  DEPTH  AND  NOMINAL 
498  m  RECEIVER  DEPTH 
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receiver  depth  at  this  site  is  in  error  and  that  the  bubble  pulse  delays 
can  be  identified  in  the  data.  To  identify  the  first  bubble  pulse,  for 
instance,  it  is  seen  that  the  negative  lobe  of  the  signal  at  25°  aligned  at 
~33  msec  is  duplicated  at  160  msec.  Also  the  arrival  aligned  at  -33  msec 
for  the  12°  signal  is  duplicated  at  -158  msec.  Thus,  the  first  bubble 
pulse  is  identified  at  -158  msec,  a  delay  of  ~125  msec,  which  is  close  to 
the  theoretical  value.  This  indicates  that  the  source  was,  in  fact, 
detonated  at  a  depth  near  18  m.  Next  it  will  be  shown  that  the  receiver 
depth  at  this  site  is  in  error.  The  source  ghost  delay  for  an  18  m  source 
is  shown  on  Fig.  36  and  the  nominal  receiver  ghost  delay  is  also  shown. 
Since  the  receiver  ghost  has  reflected  from  the  air-water  interface  one 
time,  it  is  of  a  phase  opposite  to  the  arrival  that  has  not  been  reflected 
from  the  air-water  interface.  Therefore,  since  the  first  bubble  pulse  will 
have  approximately  the  same  amplitude  and  the  same  phase  as  the  shock 
pulse,  then,  for  the  signal  at  12°,  the  arrival  at  -150  msec  should  be 
almost  canceled  to  zero.  But  it  is  not.  Also  the  absence  of  an  arrival  at 
about  230  msec  for  the  signal  at  19°  is  clear  evidence  that  the  nominal 
receiver  depth  is  in  error. 

The  arrivals  shown  in  the  box  in  Fig.  36  are  also  unexplained 
U3ing  the  nominal  source  and  receiver  depths.  These  arrivals  are  nearly 
the  same  amplitude  as  the  aligned  arrivals  and  careful  observation  leads  to 
the  conclusion  that  they  are  opposite  in  phase.  It  can  also  be  seen  for 
the  signal  at  19°  that  the  aligned  arrival  is  duplicated  with  opposite 
phase  at  about  140  msec.  Thus,  it  is  concluded  that  the  line  through  these 
arrivals  is  approximately  the  line  of  the  receiver  ghost.  This  corresponds 
to  a  receiver  depth  of  about  255  m. 
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Using  this  new,  empirically  derived  receiver  depth,  travel  time 
curves  for  the  data  with  244  m  sources  are  shown  in  Fig.  37.  This  shows 
good  agreement  with  the  data  except  for  the  early  arrival  mentioned  pre¬ 
viously.  If  the  model  of  the  ocean  bottom  as  a  specular  reflection  and  a 
refracted  component  is  accurate,  this  arrival  indicates  that  the  sound 
speed  versus  depth  increases  more  rapidly  than  the  first  order  prediction 
of  g  =  1.1  sec  -1 . 

In  summary,  this  section  has  shown  that  trends  of  acoustic  energy 
can  be  identified  (n)  travel  time  versus  grazing  angle  plots  of  acoustic 
signals  when  a  particular  phase  is  aligned.  The  interpretation  of  these 
trends  is  facilitated  by  constructing  travel  time  curves  of  the  filter 
components  of  the  signal  and  overlaying  these  travel  time  curves  on  the 
data. 

Overlaying  the  travel  time  curves  on  the  data  can  also  reveal 
inconsistencies  between  nominal  source  and/or  receiver  depths  and  a 
priori  estimates  of  geoacoustic  structure.  The  true  depth  of  the  source 
and/or  receiver  can  be  measured  from  the  aligned  data.  These 
measurements  can  be  greatly  facilitated  by  the  proper  choice  of 
experimental  geometry.  And,  in  fact,  clear  and  unambiguous 
identification  of  filter  components  is  possible  for  some  experimental 
geometries.  Good  estimates  of  geoacoustic  structure  from  aligned  data 
can  also  be  greatly  facilitated  by  carefully  chosen  experimental 
geometry. 

In  the  next  section  correlation  techniques  will  be  used  to 
measure  delays  in  the  acoustic  signals  as  a  function  of  grazing  angle. 
The  main  advantage  of  these  techniques  is  that  the  correlations  are 
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self-aligned  with  respect  to  the  zero  lag  value  and  the  delays  between 
convolutional  components  are  more  readily  seen  and  measured. 


CHAPTER  3 

CORRELATION  TECHNIQUES 


3. 1  Correlation  for  Time  Delay  Detection 

Originally,  when  this  research  was  begun,  it  was  thought  that 
the  ocean  bottom  transfer  function  could  be  described  by  the  two- 
component  model  of  a  specularly  reflected  and  a  refracted  arrival.  In 
this  case,  if  the  time  delay  between  these  components  could  be  found  as 
a  function  of  grazing  angle,  then  velocity  as  a  function  of  depth  into 
the  ocean  bottom  could  be  determined.  Since  the  power  cepstrum  had  been 
found  to  be  useful  for  time  delay  detection  of  bubble  pulses  (Mitchell 
et  al.,  1976),  it  was  thought  that  for  some  experimental  geometries  and 
geoacoustic  structures  it  would  also  be  useful  in  extracting  the  delay 
between  the  reflected  and  refracted  components. 

Later  it  was  also  decided  to  attempt  to  use  the 
autocorrelation  for  detection  of  the  impulse  response  delays.  The 
autocorrelation  has  the  advantage  that  it  has  the  superposition  property 
with  respect  to  convolution  while  the  cepstrum  does  not  fi.e., 
AC(x*y)=AC(x)*AC(y)  .  The  power  cepstrum  does  not  have  the 
superposition  property  with  respect  to  convolution  due  to  a  logarithm 
operation  (see  Appendix  II).  The  effect  of  the  logarithm  operation  in 
the  calculation  of  the  pov-;1"  cepstrum  is  to  give  peaks  at  all  sum  and 
difference  delays  as  well  as  at  the  delays  in  the  signal.  In  addition, 
there  will  be  peaks  on  the  power  cepstrum  at  all  multiples  of  the  delays 
in  the  signal  and  at  all  multiples  of  the  sura  and  difference  delays. 
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This  produces  a  complicated  pattern  of  cepstral  components  with 
amplitudes  that  are  not  readily  predicted.  Thus,  the  autocorrelation  is 
initially  preferred  over  the  power  cepstrum,  since  the  superpostion 
property  allows  for  the  amplitude  and  time  arrival  of  interfering 
components  to  be  more  easily  predicted. 

3.2  Theory  of  the  Autocorrelation 

The  autocorrelation  is  simply  an  autoconvolution  in  which 
neither  of  the  component  factors  are  reversed  before  integration. 
Figure  38  illustrates  the  value  of  the  autocorrelation  of  f(u)  at  the  point 
x.  The  autocorrelation  is  defined  by  the  integral 


Pff(x) 


f(u)  f(u-x)du 


9 


which  is  clearly  the  area  of  the  produce  of  the  function  f  with  a  shifted 
version  of  itself  by  the  amount  x.  Substituting  w=u-x,  we  have  also 


f(w+x)  f(w)dw 


so 


pff(x)  =x* 

=  f  f(u)  f(u+x)du  =  f  f(u)  f(u-x)du 

•'-Qo  •/-o • 


Thus,  it  does  not  matter  in  which  direction  the  function  is  shifted  to 
calculate  the  autocorrelation.  This  means  that,  if  f  is  real,  then  its 
autocorrelation  P^  is  even. 

Another  property  of  the  autocorrelation  function  is  that  phase 
information  is  lost.  This  is  most  easily  shown  by  considering  its  Fourier 
transform.  A  time  signal  is  a  real  function,  and  therefore  its 
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autocorrelation  is  even  and  its  Fourier  transform  degenerates  to  the 
cosine  transform  which  is  real  and  even,  i.e., 

Pff(x)  =  2^  Pff(x)  cos(27rxs)dx 
"  2il'  Pff(x)  sin(2jrxs)dx 

The  second  integral  goes  to  zero  so  the  Fourier  transform  of  the  auto¬ 
correlation  has  no  imaginary  part  and  the  phase  information  of  the  original 
function,  f,  is  lost. 

Finally,  two  other  comments  are  warranted  about  the  use  of  the 
autocorrelation  in  this  work.  First,  the  autocorrelation  has  the 
superposition  property  with  respect  to  convolution:  the  autocorrelation 
of  the  convolution  of  two  functions  is  equal  to  the  convolution  of  the 
autocorrelations  of  the  functions  (see  Appendix  II).  This  mean  that  the 
effects  of  coherent  noise  such  as  bubble  pulses  and  ghosts  are  more  easily 
recognized  and  predicted.  Second,  for  a  time  series  with  echoes  or  dis¬ 
torted  echoes,  it  would  be  expected  that  a  contribution  to  the  autocorrela¬ 
tion  will  occur  at  the  time  shifts  corresponding  to  the  echo  delays,  and, 
although  the  phase  information  is  lost,  the  effect  of  distortion  such  as  90 
or  180°  phase  shifts  are  apparent  in  the  sign  and  shape  of  the  autocorrela¬ 
tion  peaks  at  the  corresponding  time  delays. 

3«3  Theory  of  the  Power  Cepstrum 

The  power  cepstrum  is  the  power  transform  of  the  log  power 
spectrum  of  a  function.  The  significance  of  this  function  is  shown  in  the 
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following.  Consider  a  time  function  which  has  a  single  echo  (Kemerait, 
1972). 

x(t)  =  s(t)  +  a0s(t-tQ)  .  (3.3.1) 

The  power  spectrum  of  x(t)  can  be  written 

<*>x(u>)  =  0g(w)jj  +  a^  +  2aQ  cosuj  tQj  ,  (3.3.2) 

and  the  logarithm  of  both  sides  is 

log(m  (w))  =  log  0  (oj)  +  log |1  +  a^  +  2a  cosu  t  1  .  (3-3.3) 

x  s  l  o  o  oj 

For  aQ<1  the  second  term  can  be  approximated  by  2aQcosujto.  Thus,  if  s(t) 

is  an  impulse,  then  the  log  power  spectrum  is  essentially  a  cosinusoidal 

ripple  for  a  «1.  The  power  transform  of  log(4>x(u>) )  produces  a  spike  at 

the  quefrency,  tQ,  of  the  ripple  (tQ  is  the  reciprocal  period  of  the 

sinuosoid  appearing  on  the  log  power  spectrum,  measured  in  quefrencies) . 

In  the  case  where  a„  is  not  much  less  than  one,  additional  terms  of  the 

series  expansion  are  required  which  add  spikes  to  the  power  cepstrum  at 

integral  multiples  of  tQ.  Essentially,  the  same  results  apply  for  aQ»1 

2 

which  can  be  shown  by  dividing  the  second  term  of  Eq.  3.3.2  by  aQ.  This 
yields 


Vw)  =  ao  V->[,  +  4" +  ~r0  coswtl  . 

“o 

As  aQ  approaches  one,  the  presence  of  multiple  quefrencies  is  more 
apparent. 
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To  use  the  power  cepstrum  for  detecting  time  delays  within  the 
acoustic  signals  discussed  in  this  thesis,  at  least  two  complicating 
effects  must  be  considered.  First,  for  the  time  delay  of  interest  the 
delayed  arrival  has  been  distorted  by  a  90°  phase  shift.  Figure  39  shows 
that  the  effect  of  this  is  to  produce  an  even,  negative  sinusoidal  ripple 
in  the  log  power  spectrum  rather  than  a  positive  cosinusoidal  ripple. 
Since  the  power  transform  of  a  sinuosoid  simply  gives  a  spike  at  the 
frequency  (or  quefrency)  of  the  sinusoid,  the  power  cepstrum  will  be 
unaffected  by  the  90°  phase  shift  of  the  refracted  component. 

Second,  the  acoustic  signals  will  have  multiple  echoes  (i.e., 
multiple  bottom  bounce  impulse  response  components,  bubble  pulses,  and 
ghosts).  This  will  have  the  effect  of  ripples  in  the  log  power  spectrum  at 
the  quefrencies  associated  with  each  delay  and  at  all  difference  delays. 
This  is  complicated  by  the  modulating  effect  of  beating  between  these  que¬ 
frencies.  This  beating  between  quefrencies  is  the  nonlinear  effect  of  the 
power  cepstrum  operation,  giving  values  on  the  power  cepstrum  at  all  sum 
and  difference  delays. 

Although,  theoretically,  the  fact  that  the  absolute  value  of  the 
ghost  amplitudes  is  nearly  one  would  seem  to  negate  the  usefulness  of  time 
delay  detection  by  the  power  cepstrum  technique,  this  does  not  appear  to  be 
the  case  as  shown  for  the  synthetic  and  real  acoustic  signals  that  follow. 
3.4  Synthetic  Results 

In  this  section  it  is  attempted  to  model  acoustic  signals  and 
their  correlations  for  a  range  of  relative  source-receiver  geometries. 
Figure  40  is  a  flowchart  illustrating  the  format  of  the  plots  that  illu¬ 
strate  the  results.  The  acoustic  signal  is  modeled  in  time  with  the 


x(t)  =  s(t)  +  [aQ  s(t-tQ)  «  -1] 
x(t|  =  s(t|  •  (6(t|  +  aQ  5<t-tQ)  ♦ 
X(f)  =  S(f)  [1  +  (aQ  ®iwto)  (isgnf)] 


$xx<f>  =  ^ss**1  [1  +  (a0  eitt,t°>  ( isgnf ) I  [1-(a„  e'iw,o)  (isgnf)l 
=  $ss(f)  [1  +  (isgnf) (2a0)(i  sin  wtQ)  -  a*  (isgnf)2] 

=  4>$s(f)  (1  -  2aQ(sgnf)  sin  cot0  +  a*] 

!°g  $xx(f)  **  log  $55  ~  2a0  (sgnf)  sin  cotQ 
if  s(t)  =  6(t)  and  0  <  aQ  #  1 


FIGURE  39 

DERIVATION  OF  LOG  POWER  SPECTRUM  OF 
DELAYED  PULSE  WITH  90°  PHASE  SHIFT 
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TIME  SIGNATURE 


OUTPUT  PLOT 


OUTPUT  PLOT 


OUTPUT  PLOT 


FIGURE  40 

BLOCK  DIAGRAM  SHOWING  FLOWCHART  OF  THE  CALCULATION 
OF  THE  AUTOCORRELATION,  LOG  POWER  SPECTRUM,  AND 
POWER  CEPSTRUM  OF  A  SIGNAL 
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appropriate  bubble  pulses  and  ghosting  effects.  The  synthesized  waveform 
was  derived  by  the  method  discussed  in  Section  1.5. 


Figure  41  shows  the  impulse  response  at  26°  predicted  using  the 
geoacoustic  model  shown  in  Table  II.  It  is  seen  that  there  are  two  main 
components  and  a  smaller  component.  The  specular  reflection  is  the  first 
arrival,  with  an  amplitude  of  about  0.25  of  the  incident  wave.  The 
refracted  component  has  a  90°  phase  shift  and  has  undergone  some  attenua¬ 
tion.  The  third  component  is  a  refracted  wave  which  has  reflected  once 
from  the  underside  of  the  water-sediment  interface,  has  undergone  a  total 
of  a  360°  phase  shift,  and  has  undergone  twice  the  attenuation  of  the 
refracted  arrival  (recall  Fig.  16).  The  specular  reflection  was  aligned 
at  25  msec. 

The  autocorrelation  shown  beneath  the  time  signal  reveals  that 
an  arrival  which  is  similar  in  shape  to  the  refracted  waveform  occurs  in 
the  autocorrelation.  The  middle-zero  crossing  of  this  arrival  is  the  time 
delay  between  the  reflected  and  refracted  components.  Note  that  the 
zero  lag  value  of  the  autocorrelation  is  aligned  at  25  msec. 

The  log  power  spectrum  beneath  the  autocorrelation  has  some 
interesting  features.  First,  the  dropoff  of  power  with  increasing  fre¬ 
quency  is  essentially  the  dropoff  due  to  attenuation  along  the  refracted 
path.  It  is  seen  that  there  is  also  a  very  low  quefrency  ripple  on  the 
log  power  spectrum  which  overlies  the  linear  dropoff  due  to  attenuation. 
This  low  quefrency  component  may  be  due  to  the  apparent  time  delay 
between  the  positive  and  negative  lobes  of  the  refracted  arrival. 

The  amplitude  of  the  higher  quefrency  ripple  in  the  log  power 
spectrum  is  frequency  dependent.  This  is  not  unexpected  because 
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TABLE  II 


GEO ACOUSTIC  MODEL  FOR  FIGS.  41-48 


Depth 

Velocity 

(m) 

(tn/sec ) 

o- 

1515 

0+ 

1515 

1000“ 

2515 

1000+ 

5200 

Attenuation 

Density 

(dB/m-kHz) 

(K/cm3) 

1.0 

0.02 

1.7 

0.02 

2.4 

0.01 

2.5 
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attenuation  has  been  included  in  the  model.  At  low  frequencies  the 
refracted  energy  is  not  attenuated  and  aQ  is  much  greater  than  one  and  the 
ripple  is  small.  At  higher  frequencies  the  attenuation  is  greater,  aQ 
approaches  one,  and  the  amplitude  of  the  ripple  increases.  Eventually,  a 
frequency  is  reached  where  the  refracted  energy  is  attenuated  below  the 
level  of  the  reflected  arrival.  At  this  point  the  amplitude  of  the  ripple 
begins  to  decrease  since  it  is  directly  proportional  to  aQ,  which 
decreases. 

The  power  cepstrum  of  the  signal  shown  in  Fig.  41  shows  that  the 
time  delay  of  interest  is  easily  detected.  Zero  quefrency  is  located  at 
25  msec.  The  value  of  the  time  delay  will  be  at  the  center  of  the  two  peaks 
at  about  90  msec  (tQ  =  65  msec).  Large  cepstral  peaks  would  also  occur  at 
the  low  quefrencies  due  to  attenuation,  but  have  been  muted  for  the  pur¬ 
poses  of  detecting  the  reflection-refraction  d^ay.  The  first  three  sam¬ 
ples  of  the  power  cepstrum  are  zeroed  on  this  and  all  other  power  cepstrum 
plots  with  the  exception  of  Fig.  -53- 

Figure  42  shows  the  effect  of  adding  a  244  m  source  to  the 
previous  impulse  response.  It  is  seen  that,  although  the  impulse  response 
is  highly  distorted,  the  time  delay  of  interest  at  about  90  msec  in  the 
autocorrelation  occurs  with  little  distortion. 

The  effect  of  the  source  on  the  log  power  spectrum  and  power 
cepstrum  is  shown  in  Fig.  42  as  well:  significant  new  ripples  associated 
with  the  first  and  second  bubble  pulses  are  shown  in  the  log  power  spec¬ 
trum.  These  ripples  show  up  as  spikes  on  the  power  cepstrum  plot  at  the 
bubble  pulse  delays  of  about  19  and  33  msec.  No  significant  beating 
between  quefrencies  is  observed,  but  the  amplitude  of  the  peak  associated 


80 


TIME  SIGNATURE 


AUTOCORRELATION 


i - ' - 1 - - i — - — i - i - - — i 

0  100  200  300 

TIME  —  msec 


POWER  CEPSTRUM 


0 


WO  ~ 1  200 

QUEFRENCY  -  msec 


1 - 1 

300 


FIGURE  42 

TIME  SIGNATURE,  AUTOCORRELATION,  POWER  SPECTRUM,  AND 
POWER  CEPSTRUM  OF  PLANE  WAVE  IMPULSE  RESPONSE 
CONVOLVED  WITH  SOURCE  WAVEFORM  (244  m  SOURCE) 


81 


ARL:UT 
AS-8 1-929 
GDI  GA 
8  5-81 


with  the  refracted  arrival  is  diminished  from  the  simple  impulse  response 
for  the  power  cepstrum. 

Figure  43  shows  the  effect  of  adding  the  244  m  ghost.  No 
significant  distortion  is  added  to  either  the  time  or  autocorrelation 
plots.  The  time  delay  is  equally  well  defined  in  this  figure  as  compared 
to  Fig.  42.  On  the  other  hand,  the  high  quefrency  ripple  associated  with 
the  source  ghost  dominates  both  the  log  spectrum  and  power  cepstrum.  The 
power  cepstrum  of  this  figure  indicates  that  trace  to  trace  adjustment  of 
amplitude  may  be  necessary  to  maximize  the  visibility  of  the  desired  delay. 
Unfortunately,  the  optimum  trace  to  trace  adjustment  is  probably  difficult 
to  predict  since  it  is  due  to  the  complicated  beating  pattern  of  que- 
frencies  associated  with  the  experiment  geometry  and  geoacoustic  structure 
of  the  ocean  bottom.  In  addition  to  the  reduction  of  amplitude  of  the 
desired  delay  peak  due  to  the  beating  of  quefrencies,  there  appears  to  be  a 
relative  amplification  of  cepstral  peaks  other  than  the  reflected- 
refracted  peak. 

Figure  44  shows  the  effect  of  including  a  receiver  ghost  and  a 
source  plus  receiver  ghost  to  the  signal.  In  this  case  a  shallow  (75  m) 
receiver  was  chosen  to  observe  the  effect  of  adding  a  ghost  between  the 
reflected  and  refracted  components.  It  is  seen  that  this  geometry  (or 
arrival  structure)  considerably  obstructs  identif ication  of  the  reflected- 
refracted  delay.  For  the  time,  autocorrelation,  and  log  spectrum  dis¬ 
plays,  it  is  virtually  impossible  to  detect  the  reflected-refracted  delay- 
For  the  power  cepstrum  this  is  very  unlikely,  but  perhaps  on  the  two- 
dimensional  T -6  plot  the  arrival  could  be  identified.  The  two  boxed 
arrivals  are  of  particular  concern  because  they  are  nearly  as  large  as  the 
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refracted  arrival,  but  have  no  obvious  explanation  in  terms  of  signal 
filter  components.  Actually  they  are  beats  between  the  other  quefrency 
components . 

Since  we  often  study  acoustic  signals  that  have  interacted  with 
the  bottom  a  number  of  times,  Fig.  45  illustrates  the  effect  of  a  five 
bottom  bounce  impulse  response.  There  is  no  ambiguity  of  the  time  delay 
associated  with  the  reflected  and  refracted  components  in  any  of  the  fig¬ 
ures.  The  greater  attenuation  due  to  more  bottom  refractions  is  evident  on 
the  time,  autocorrelation,  and  log  power  spectrum.  Note  that  the  time 
signature  and  autocorrelation  plots  are  smeared  in  time  and  the  dynamic 
range  of  the  log  power  spectrum  is  increased  by  five  times.  The  attenua¬ 
tion  due  to  the  five  bottom  interactions  of  the  signal  shown  in  Fig.  45 
provides  a  good  illustration  of  the  major  advantage  of  the  power  cepstruro 
technique  over  the  autocorrelation.  This  advantage  is  that  the  power 
cepstrum  gives  a  spike  at  the  reflection-refraction  delay  rather  than  a 
pulse  with  the  crosspower  spectrum  of  the  reflected  and  refracted  compo¬ 
nents.  Thus,  the  power  cepstrum  has  much  greater  resolving  power  of 
individual  impulse  response  components  than  the  autorcorrelation. 

Figure  46  shows  the  effect  of  a  five  bottom  bounce  signal  with 
the  source  and  receiver  depths  as  in  Fig.  44.  *  seen  that  the  capa¬ 
bility  of  measuring  the  desired  delay  is  reduced  U.  the  time  plot,  but 
dramatically  improved  for  the  autocorrelation  and  power  cepstrum.  The  low 
pass  earth  filtering  apparently  has  bolstered  the  reflection-refraotion 
delay  at  the  expense  of  the  bubble  pulse  delay  in  the  autocorrelation. 
This  is  due  to  greater  high  frequency  content  in  the  bubble  pulse  than  in 
the  refracted  arrival.  So,  since  attenuation  is  linearly  dependent  on 
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frequency,  low  frequencies  are  preferentially  passed  for  each  bottom 
interaction  while  the  high  frequencies  are  attenuated. 


The  amplitude  of  the  reflection-refraction  delay  is  also  greatly 
amplified  with  respect  to  other  components  on  the  power  cepstrum  plot  as 
well.  This  is  due  to  the  fact  that  for  synthetic  signals  we  have  infinite 
dynamic  range,  and  adjacent  frequency  bands  are  affected  by  the  attenua¬ 
tion  in  the  subbottom  five  times  as  much  for  the  five  bounce  signal  as  the 

C 

one  bounce  signal  (log  R  (w)  =  5  log  R(o>)).  Thus,  the  amplitude  of  the 
ripples  associated  with  the  refleeted-refracted  delay  is  five  times 
greater  for  the  five  bottom  bounce  signal  than  the  one  bottom  bounce 
signal,  and  this  effect  is  seen  in  the  power  cepstrum. 

Figure  47  shows  the  effects  on  the  time  plots  and  frequency  plot 
of  a  five  bottom  bounce  signal  with  a  244  m  source  and  a  400  m  receiver. 
This  greatly  improves  the  ability  to  detect  the  refleeted-refracted  time 
delay  in  both  the  time  and  autocorrelation  plots.  The  ability  of  the  power 
cepstrum  to  measure  the  time  delay  of  interest  is  about  the  same.  In 
general  the  ability  of  the  autocorrelation  and  power  cepstrum  to  detect 
delays  between  impulse  response  components  is  greatly  dependent  on  the 
source-receiver  geometry. 

Finally,  Fig.  48  shows  the  time  plots  and  frequency  plot  for  the 
source-receiver  geometry  empirically  verified  in  Chapter  2  for  some  real 
data.  It  is  shown  that  the  refleeted-refracted  time  delay  is  clearly 
identifiable  in  all  plots. 

The  geoacoustic  structure  also  affects  the  ability  of  the 
autocorrelation  and  power  cepstrum  for  time  delay  detection.  Figure  49 
shows  an  impulse  response  at  26°  with  an  ocean  bottom  with  the  geoacoustic 
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structure  shown  in  Table  III  rather  than  that  of  Table  II.  The  higher 
sound  speeds  result  in  shorter  pathlengths  of  the  refracted  energy  and  the 
lower  density  at  the  water-sediment  interface  results  in  the  specular 
reflection  being  smaller.  Thus,  the  refracted  event  is  larger  relative  to 
the  specular  reflection  in  Fig.  49  compared  to  Fig.  41.  Comparing  the 
autocorrelation  and  power  cepstrum  plots  of  Figs.  41  and  49,  it  is  seen 
that  the  change  in  impulse  responses  does  not  affect  the  ability  of  the 
autocorrelations  or  power  cepstrums  in  detecting  the  time  delay  of  inter¬ 
est,  although  the  log  power  spectrum  Indicates  considerably  greater 
dynamic  range  for  the  impulse  response  of  Fig.  41.  The  effect  of  this 
difference  in  dynamic  range  on  the  power  cepstra  would  be  seen  if  they 
were  not  self-normalized.  The  power  cepstra  in  these  plots  are  normalized 
after  the  first  three  samples  and  last  three  samples  are  zeroed. 

Figure  50  shows  a  5  BB  synthetic  acoustic  signal  with  the  same 
source  and  receiver  geometry  as  the  data  discussed  in  Chapter  2  using  the 
impulse  response  shown  in  Fig.  49.  It  is  seen  that  the  autocorrelation  and 
power  cepstrum  readily  reveal  the  time  delay  of  interest. 

Comparing  Fig.  48,  which  has  the  same  source-receiver  geometry 
as  Fig.  50,  but  uses  an  impulse  response  like  that  of  Fig.  41  rather  than 
Fig.  49,  it  is  seen  that  there  are  significant  differences.  In  particular, 
the  bubble  pulse  delays  and  the  source  and  receiver  ghost  delays  of  the 
power  cepstrum  plot  in  Fig.  50  are  much  stronger  than  those  in  Fig.  48 
relative  to  the  reflection-refraction  delay.  This  is  due  to  the  consider¬ 
ably  greater  dynamic  range  of  the  reflection-refraction  quefrency  in  the 
log  power  spectrum.  The  greater  dynamic  range  of  the  log  power  spectrum  of 
Fig.  48  is  due  to  attenuation  along  the  refracted  path.  The 
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autocorrelation  of  Fig.  48  also  isolates  the  reflection-refraction  delay 
better  than  does  the  autocorrelation  of  Fig.  50.  This  is  also  due  to  the 
longer  refracted  pathlengths  in  the  geoacoustic  model  of  Fig.  48,  which 
results  in  the  power  being  more  concentrated  at  low  frequencies  where  the 
refracted  energy  is  strong  relative  to  the  bubble  pulses. 

Next  the  effect  of  white  uncorrelated  noise  will  be  examined  on 
the  power  cepstrum.  This  will  be  done  by  adding  some  small  fraction  to 
the  zero  lag  value  of  the  autocorrelation.  This  fraction  will  represent 
the  percentage  of  the  power  of  the  time  signal  which  was  made  up  of 
white  uncorrelated  noise.  This  only  contributes  to  the  autocorrelation 
integral  at  the  zero  lag  value. 

Figure  51  shows  the  effect  of  adding  0.01  to  the  zero  lag 
value  of  the  autocorrelation.  It  is  seen  that  the  effect  of  this  is  to 
reduce  the  dynamic  range  of  the  log  power  spectrum  and  this  has  a 
corresponding  effect  on  the  power  cepstrum  of  reducing  the  amplitude  of 
the  peak  due  to  the  reflection-refraction  delay  and  the  bubble  pulse 
delays.  This  is  due  to  the  fact  that  these  power  cepstrum  peaks  are 
more  dependent  on  the  dynamic  range  at  the  high  frequencies  than  the 
ghosts  are,  and  the  fact  that  the  high  frequencies  are  below  the  mean 
power  level  of  the  spectrum.  Figure  52  shows  the  effect  of  having  10  dB 
signal  power  excess.  It  is  seen  that  the  level  of  the  reflection- 
refraction  peak  and  bubble  pulse  peaks  are  further  reduced  and 
distorted. 

Finally,  in  this  section  the  common  procedure  of  whitening  the 
log  power  spectrum  before  calculating  the  power  cepstrum  is  examined  for 
synthetic  plane  wave  data.  This  is  accomplished  by  removing  the  slope 
and  mean  of  the  log  power  spectrum  by  linear  regression.  Figure  53 
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FIGURE  61 

PLOT  SHOWING  THE  EFFECT  OF  ADDING  WHITE 
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effect  of  whitening  the  log  power  spectrum  before  calculating  the  power 
cepstrum  where  the  first  three  samples  of  the  power  cepstrum  are  not 
zeroed.  It  is  seen  that  the  bubble  pulse  delay,  the  reflection-refraction 
delay,  and  the  ghost  delays  can  be  detected.  The  significance  of  this  is 
that,  ordinarily,  without  whitening  the  log  power  spectrum  the  first  three 
samples  of  the  power  cepstrum  totally  dominate  the  plot.  Therefore,  one 
alternative  to  zeroing  the  first  three  samples  of  the  power  cepstrum  is  to 
whiten  the  log  power  spectrum  before  calculating  the  cepstrum. 

Figure  54  shows  the  effect  of  whitening  the  log  power  spectrum 
before  calculating  the  power  cepstrum  and  zeroing  the  first  three  samples 
of  the  power  cepstrum.  A  comparison  of  this  figure  with  Fig.  50  indicates 
that,  for  this  synthetic  example,  whitening  the  log  power  spectrum  has 
little  effect  on  the  power  cepstrum. 

3.5  Data  Autocorrelogram  and  Power  Cepstra 

In  this  section  autocorrelations  and  power  cepstra  are 
calculated  for  real  data  obtained  in  a  turbidite  province.  Figure  55  shows 
an  autocorrelogram  with  the  same  source-receiver  geometry  as  the  data 
shown  in  Fig.  33.  The  only  differences  between  the  data  shown  in  Fig.  33 
and  Fig.  55  are  that  the  data  for  the  autocorrelogram  are  sampled  almost 
every  degree  and  extend  to  40°  rather  than  30°. 

The  source  and  receiver  ghost  delays  and  the  source  bubble  pulse 
delays  are  evident  on  the  autocorrelogram.  The  depth  of  detonation  of  the 
source  has  some  variance  and  this  is  seen  in  the  variance  of  the  bubble 
pulse  delay  and  source  ghost  delay.  For  these  data  even  the  second  bubble 
pulse  is  observable  on  the  autocorrelations,  particularly  at  low  grazing 

angles.  Since  the  receiver  is  fixed,  the  receiver  ghost  delay  is  a 
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FIGURE  64 

PLOT  SHOWING  THE  EFFECT  OF  WHITENING  THE  LOG  POWER  SPECTRUM 
BEFORE  CALCULATING  THE  POWER  CEPSTRUM  WITH  THE  FIRST  THREE 
SAMPLES  OF  THE  POWER  CEPSTRUM  ZEROED 
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AUTOCORRELOGRAM  FOR  A  TURBIDITE  PROVINCE 
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constant.  The  receiver  ghost  delay,  in  particular,  and  all  other  delays 
can  best  be  seen  by  looking  up  the  page  with  your  eye  level  to  the  page. 

The  curved  lines  drawn  on  the  autocorrelation  plot  represent  the 
time  delay  between  reflected  and  refracted  energy.  Currently,  it  is  not 
known  whether  this  time  delay  between  the  reflected  and  refracted  energy  is 
due  to  a  specular  or  subbottoa  reflection.  It  is  seen  that  this  delay  is 
well  defined  between  about  15-28°.  Beyond  about  32°  there  is 
interference  between  the  reflected-refraeted  time  delay  of  the  zero  lag 
value  and  the  ghost  refracted  time  delay.  This  ghost  refracted  delay  is 
the  delay  between  the  refracted  component  of  the  first  arrival  and  the 
reflected  arrival  of  the  ghost.  This  delay  bends  baok  toward  the 
ordinate  and  the  event  has  the  same  shape  as  the  reflected-refraeted 
delay  of  the  zero  lag  value. 

The  shape  of  the  ghost  events  themselves  are  the  negative  of 
the  zero  lag  value;  however,  this  is  not  clearly  seen  in  Fig.  55  since 
the  source  and  receiver  depths  are  nearly  the  same.  The  souroe  plus 
receiver  ghost  energy  is  also  seen  in  Fig.  55. 

At  2.5°  the  sum  of  all  of  the  autocorrelations  is  plotted. 
This  provides  an  estimate  of  the  autocorrelation  of  the  source  waveform 
assuming  that  attenuation  is  not  significant,  and  that  all  events  with 
other  than  zero  time  moveout  as  a  function  of  grazing  angle  cancel  to 
zero. 

A  significant  improvement  in  the  autooorrelogram  for  measuring 
the  refleoted-refracted  delay  is  obtained  by  muting  the  time  waveforms 
beyond  the  source  ghost  delay  of  the  first  break.  This  is  beoause  the 
bottom  angles  for  the  ghosted  arrivals  are  greater  than  the  arrival 

which  has  not  been  ghosted.  Thus,  the  impulse  response  delays  on  the 
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autocorrelations,  which  have  not  had  the  tine  signals  muted  beyond  the 
source  ghost,  give  poorer  resolution  of  the  reflected-refracted  event. 
This  increased  resolution  of  the  reflected-refracted  event  is  seen  in 
Figs.  55  and  56. 

Figure  57  shows  an  autocorrelogram  for  18  n  shot  data  with  a 
receiver  depth  of  256  m  in  the  same  turbidite  province  as  the  244  m  shot 
data.  The  source  and  receiver  ?  -'&cs  are  identifiable  and  the  bubble 
pulses  are  apparent  from  tb*  '  -ace.  The  pulse  at  about  150  msec 

is  the  crossoorrelation  of  Cfct-  pulse  with  the  first  bubble  pulse, 

the  pulse  at  about  200  is  l he  crosscorrelation  of  the  shook  pulse 

with  the  second  bubble  puX*e,  and  the  pulse  at  about  135  msec  is  the 
crosscorrelation  of  the  first;  and  second  bubble  pulse.  The  variability 
of  the  bubble  pulse  delays  indicates  the  variance  of  the  depths  of 
detonation  of  the  source. 

The  reflected-refracted  delay  of  the  18  m  source  auto- 
correlograms  is  the  same  as  that  on  the  244  m  source  autocorrelograms 
which,  of  course,  is  expected.  The  18  m  shot  data,  however,  also  reveal 
significant  reflectivity  not  resolved  by  the  244  m  data.  It  is 
difficult  to  resolve  more  than  one  reflector  since  this  is  five  bottom 
bounce  data  and  some  of  the  peaks  are  multiples,  but  a  reflector  which 
is  delayed  slightly  from  the  zero  log  value  is  well  defined.  Thus, 
there  appears  to  be  a  significant  shallow  refleotor  in  this  turbidite 
province.  To  improve  the  resolution  of  the  reflected-refracted  delay 
the  time  signals  are  muted  past  the  receiver  ghost  before 
autocorrelating  in  Fig.  58. 

Returning  to  Fig.  57,  it  is  seen  that  at  angles  of  16,  18,  and 

21°  the  main  components  of  the  autocorrelation  are  reproduced  for 
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different  shots — which  is  reassuring.  Comparing  autocorrelations  at  the 
same  angle  reveals  the  approximate  noise  level  of  a  single 
autocorrelation. 


One  method  of  reducing  the  noise  level  of  the  autocorrelation 
and  at  the  same  time  reducing  the  effect  of  bubble  pulses  is  to  sum 
autocorrelations  at  the  same  angle  for  different  numbers  of  bottom 
bounces.  An  example  of  this  is  shown  for  a  bottom  angle  of  26°  in 
Fig.  59.  It  is  seen  here  that  there  are  similarities  of  the 
autocorrelations  and  differences.  The  similarities  are  due  to  the  fact 
that  the  bottom  angle  and  source-receiver  geometries  are  nearly  the 
same.  The  dissimilarities  are  due  to  source  depth  differences,  slightly 
different  bottom  angles,  and  different  bounce  points  on  the  bottom. 
Autocorrelations  like  that  for  the  8  BB  signal  are  interpreted  to  be  due 
to  bad  bounce  points.  The  summed  autocorrelation  is  shown  at  the  bottom 
of  Fig.  59.  The  peaks  in  the  summed  autocorrelation  are  readily 
interpreted.  The  first  positive  peak  is  likely  the  bubble  pulse  at  the 
mean  source  depth,  the  second  event  is  the  reflection-refraction  delay 
with  its  bubble  at  its  tail,  the  third  event  is  some  combination  of  the 
ghost  delays  and  the  multiple  of  the  reflection-refraction  delay,  and 
the  waveform  beyond  the  third  event  is  associated  with  ghost  energy. 
Thus,  it  appears  that  an  autocorrelogram  of  summed  autocorrelations 
would  significantly  facilitate  time  delay  detection  of  events  in  the 
acoustic  signal. 

Figures  60-63  show  examples  of  the  power  cepstrum  applied  to 
the  data.  Figure  60  shows  the  power  cepstrums  of  the  data  shown  in 
Figs.  55  and  56.  The  first  three  samples  are  zeroed  for  display 
purposes  and  zero  quefrency  is  located  at  0  msec.  The  first  bubble 

pulse  is  the  strongest  event  in  Figs.  60  and  61  and  the  source  and 
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FIGURE  SHOWING  AUTOCORRELATIONS  WITH  A 
BOTTOM  ANGLE  OF  ABOUT  26°  FOR  DIFFERENT 
NUMBERS  OF  BOTTOM  BOUNCES  AND  THEIR  SUM 
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FIGURE  60 

POWER  CEPSTRUM  OF  DATA  IN  TURBIDITE  PROVINCE 
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FIGURE  63 

SIGNAL  MUTED  PAST  SOURCE  GHOST  DELAY 
AND  POWER  SPECTRUM  WHITENED  BEFORE 
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receiver  ghost  energy  is  also  evident.  The  source  plus  receiver  ghost 
energy  is  less  evident.  The  fact  that  the  power  oepstrum  does  not  have 
the  superposition  property  with  respect  to  convolution  accounts  for  the 
large  amplitude  variability  of  the  ghost  arrivals.  This  variability  is 
due  to  beating  between  all  quefrencies  of  the  signal.  The  quefrencies 
that  are  most  greatly  affected  will  be  determined  by  the  amplitudes  and 
periodicity  of  the  quefrencies  which  occur.  Thus,  interpretation  of  a 
power  cepstrum  or  suite  of  power  cepstra  will  be  less  exact  than  the 
autocorrelation,  and  will  probably  require  some  experience.  Coherent 
noise  will  be  less  identifiable  and  the  shape  of  the  cepstral  peaks  will 
be  less  interpretable  in  terms  of  the  phase  of  the  time  arrivals  than 
for  the  autocorrelation. 

Figure  61  shows  the  power  cepstra  shown  in  Fig.  60  where  the 
power  spectra  were  whitened  before  the  power  cepstra  were  calculated. 
Some  improvement  seems  to  have  occurred.  The  reflection-refraction 
delay  is  as  easy  to  detect  and  some  subbottom  reflectivity  is  also 
apparent  at  delays  less  than  the  reflection-refraction  delay. 

Figures  62  and  63  show  the  effect  of  muting  the  time  signals 
past  the  source  ghost  delay  before  calculating  the  power  cepstra.  The 
reflection-refraction  delay  is  easier  to  detect  particularly  when  the 
power  spectra  are  whitened  before  calculating  the  power  cepstra  in 
Fig.  63.  In  these  figures  the  number  of  cepstral  peaks  at  delays  less 
than  25  msec  is  indicative  of  subbottom  reflectivity;  however, 
identification  of  events  is  not  certain.  The  inability  to  identify 
reflection  events  in  these  plots  may  be  due  to  beating  between  quefrency 
components,  or  to  an  increase  in  the  variance  of  the  noise  speotrum 

caused  by  muting  beyond  the  source  ghost  delay  of  the  first  break. 
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CHAPTER  4 


IMPULSE  RESPONSE  ESTIMATION  USING  MULTIPLE  BOTTOM  BOUNCE  SIGNALS 


4.1  Impulse  Response  Estimation 

An  alternative  approach  to  correlation  for  extracting  geoacoustic 
information  about  the  ocean  bottom  from  acoustic  data  is  the  deconvolution 
of  unwanted  signal  components  (i.e.,  estimation  of  the  ocean  bottom  impulse 
response  by  removing  bubble  pulses  and  ghosts).  Debubbling  has  been  accom¬ 
plished  for  over  ten  years  by  frequency  domain  deconvolution.  Hovem  first 
introduced  deconvolution  of  bubble  pulses  by  measuring  the  source  waveform 


and  dividing  the  source  spectrum  with  the  spectrum  of  the  bottom  interacting 
signal  (Hovem,  1979;  Hastrup,  1969).  This  is  described  by  the  equation 


I(w)  B  f{S)  x  G(w) 


(4.1.1) 


where 

I(w)  =  impulse  response  estimate, 

A(w)  =  bottom  interacting  signal, 

S(w)  =  source  waveform  spectrum,  and 
G(«*>)  s  shaping  filter  taking  band  limitation  and  noise 
into  account. 

Santaniello  et  al.  improved  this  technique  by  multiplying  both  the 
numerator  and  denominator  of  (4.1.1)  by  the  complex  conjugate  of  the  source 
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spectrum,  S*(w).  This  is  equivalent  to  match  filtering,  which  reduces  the 
effect  of  uncorrelated  noise.  In  addition,  Santaniello  et  al.  (1979)  added 
white  noise  to  the  source  power  spectrum  to  stabilize  the  transfer  function 
by  preventing  denominator  values  being  near  zero, 

k„)  .  msiM.  *  0(„)  .  (4.i2) 
|s(">)|Z.C 

Dicus  (1981)  has  improved  on  this  technique  by  minimizing  the  squared  error 
of  the  impulse  response  estimate  for  each  frequency  component.  Dicus  also 
shows  that  the  deconvolution  can  be  accomplished  using  a  single  measured 
source  waveform  by  time  scaling  to  account  for  source  depth  differences 
between  shots. 

Time  domain  approaches  to  deconvolution  have  been  in  existence  in 
the  oil  industry  for  many  years.  These  techniques  generally  employ  Weiner 
filtering  that  solves  for  filter  coefficients  in  the  least  mean  square  error 
sense  which  relates  input  to  a  desired  output.  This  generally  involves 
shaping  a  seismic  field  wavelet  to  an  interpretive  wavelet  where  the  field 
wavelet  is  estimated  statistically  (Backus  et  al . ,  1977).  The  filter 
coefficients  can  be  obtained  using  the  Levinson  algorithm  and  the  decon¬ 
volution  can  be  accomplished  either  in  the  time  or  frequency  domain  (Lackoff 
and  LeBlanc,  1975).  Wood  et  al.  (1978)  show  that  debubbling  can  also  be 
accomplished  using  the  Weiner  approach  in  which  an  inverse  filter  is  derived 
from  a  measured  source  waveform.  The  technique  which  follows  is  a  Weiner 
approach  to  obtaining  an  estimate  of  the  ocean  bottom  impulse  response 
utilizing  multiple  bottom  bounce  waveforms  that  have  equal  source  and 
receiver  depths  and  equal  bottom  grazing  angles.  This  deconvolution 
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technique  will  be  termed  multiple  bottom  bounce  deconvolution  (MBB 
deconvolution) . 

4.2  Theory  of  MBB  Deconvolution 

The  method  used  in  this  work  to  recover  an  estimate  of  the  impulse 
response  of  the  ocean  bottom  is  essentially  the  same  as  the  debubbling 
technique  employed  by  Wood  et  al.  (1978)  except  that  the  source  waveform  is 
replaced  by  an  (n-1)  bottom  bounce  waveform  with  the  same  bottom  angle  and 
source  and  receiver  depths  as  an  (n)  bottom  bounce  waveform. 

The  flowchart  reproduced  from  Wood  et  al.  in  Fig.  64  shows  the 
derivation  of  the  deconvolution  filter  D(Z).  This  filter  consists  of  the 
crosscorrelation  of  the  (n-1)  bottom  bounce  signal  with  the  autocorrelation 
of  the  zero  delay  Weiner  inverse.  Clearly  one  would  obtain  a  zero  delay 
impulse  if  the  deconvolution  operator  was  crosscorrelated  with  the  (n-1 ) 
bottom  bounce  received  waveform. 

D(Z)Sn-1(Z)  =  B(Z)  B(Z"1)  Sn_1(Z)  Sn_1(Z_1)  =  1 

There  is  no  restriction  on  the  phase  of  Sn_1(Z)  and  ghosts  may  be  included. 

This  deconvolution  technique  can  be  extended  to  derive  an 
estimate  of  the  ocean  bottom  impulse  response  by  deconvolving  an  (n-1) 
bottom  bounce  waveform,  from  an  (n)  bottom  bounce  waveform  provided  both 
waveforms  have  equal  source  and  receiver  depths  and  equal  bottom  grazing 
angles.  This  estimate  involves  an  assumption  that  the  ghost  delays  are 
dependent  only  on  the  source  and  receiver  depths  and  the  bottom  grazing 
angles.  This  becomes  more  nearly  true  as  the  number  of  bot' om  bounces 
beoomes  large. 
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FIGURE  64 

FLOWCHART  FOR  CALCULATION  OF  DECONVOLUTION  FILTER,  D{Z) 
X(Z)  IS  THE  (n  -1)  BOTTOM  BOUNCE  WAVEFORM 

[Reproduced  from  Wood  et  el.,  (1978)] 
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Assuming  that  the  conditions  above  are  met  and  that  the  ghost 
delays  are  nearly  identical  an  estimate  of  the  impulse  response  can  be 
calculated  by  computing  an  inverse  to  the  (n-1)  bottom  bounce  signal  in  the 
least  mean  square  error  sense.  This  can  be  accomplished  by  crosscorrelating 
the  (n-1)  bottom  bounce  signal  with  the  autocorrelation  of  its  zero  lag 
Weiner  inverse  and  crosscorrelation  of  this  filter  with  the  (n)  bottom 
bounce  signal. 

This  can  be  written  in  the  Z-domain 

r  .  (z)  =  sn_1(z_1)  sn(z)  =  sH-’cz"1)  S^CZ)  I(Z)  , 

n-l  ,n 

R  .  (Z)  D(Z)  =  B(Z)  B(Z~')  S0-1^1)  Sn_1(Z)  I(Z) 

n-1  ,n 

Since  B(Z)  B(Z_1)  s"'1^1)  Sn_1(Z)  »  1  then  R„  ,  (Z)  D(Z)  *  I(Z),  where 

n- 1  f  n 

I(Z)  is  the  first  bottom  bounce  impulse  response  of  the  (n)  bottom  bounce 
signal. 

To  verify  that  I(Z)  is  the  impulse  response  of  the  first  bottom 
bounce  of  the  (n)bottom  bounce  signal,  recall  the  filter  representation  of 
the  (n)bottom  bounce  acoustic  signal 

Sn(f)  =  S»(f)  Gg(f)  GR( f )  In(f) 

where 

S*(f)  =  frequency  spectrum  of  source 

G  (f)  =  source  ghost  filter 

s 

Op(f)  a  receiver  ghost  filter 
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In(f)  =  (n)  bottom  bounce  impulse  response. 

Tf  the  ghost  delays  are  assumed  to  be  identical* 

Sn~1(f)  =  S*(f)G  (f)G_(f)In_1(f) 

S  r 

Clearly  if  the  inverse  function  to  Sn”1(f)  is  applied  to  Sn(f)  the  result 
will  be  an  estimate  of  1(f). 

If  the  experiment  geometry  is  set  up  as  shown  in  Fig.  65  the 
deconvolved  trace  is  the  first  bottom  bounce  impulse  response  of  the  (n) 
bottom  bounce  signal  and  thus  the  subbottom  characteristics  could  be  mapped 
in  space.  If  the  impulse  response  can  be  assumed  to  be  space  stationary, 
ensemble  averaging  can  be  used  to  improve  the  estimate  of  the  impulse 
response  versus  grazing  angle. 

4.3  Examples  of  MBB  Deconvolution  Technique 

Figures  66  and  67  show  examples  of  the  deconvolution  technique. 
In  Fig.  66  an  (n-1)  bottom  bounce  synthetic  waveform  is  deconvolved  from  an 
(n)  bottom  bounce  synthetic  waveform  with  equal  bottom  grazing  angles  and 
equal  ghost  delays.  It  is  seen  that  the  impulse  response  is  almost  per¬ 
fectly  recovered.  The  fact  that  it  was  not  perfectly  recovered  is  because 
the  operator  length  was  limited  to  the  length  used  in  real  data.  In 
practice  an  operator  length  is  generally  chosen  to  be  between  100  and 
250  msec.  In  the  absence  of  noise  an  infinite  operator  length  would  be 
desirable  due  to  the  ghosts  (which  have  infinite  inverses);  however,  in  the 
presence  of  noise  there  is  a  tradeoff  between  deconvolving  predictable 
energy  and  loss  of  signal-to-noise  ratio  (S/N).  The  method  of  predictive 
deconvolution  (Peacock  and  Treitel,  1969)  is  well  suited  for  some  of  the 
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FIGURE  65 

GEOMETRY  FOR  DECONVOLVING  AN  (n-1)  BOTTOM  BOUNCE  WAVEFORM 
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FIGURE  66 

RECOVERY  OF  THE  IMPULSE  RESPONSE 
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FIGURE  67 
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geometries  of  interest;  however,  the  method  described  here  is  preferable 
because  of  its  ease  of  implementation. 


Figure  67  is  an  example  application  for  real  data.  It  is  seen 
that  there  is  a  considerably  greater  noise  level  in  the  deconvolved  trace 
for  real  data.  This  is  primarily  due  to  the  breakdown  of  the  assumption 
that  the  sources  were  detonated  at  the  same  depth;  but  perhaps  it  is  also 
due  to  the  breakdown  of  the  assumptions  that  the  ghost  delays  are  identical 
and  that  the  angles  of  the  (n)  bottom  bounce  arrivals  are  identical  to  the 
angles  of  the  (n-1)  bottom  bounce  arivals.  In  spite  of  the  fact  that  these 
assumptions  will  always  break  down  to  some  extent,  Fig.  67  shows  an  all-pass 
example  in  which  the  refracted  arrival  is  identifiable  and  at  least  two 
reflectors  are  apparent.  Behind  the  refracted  arrival  there  appears  to  be 
another  event  which  may  be  either  a  deeper  refraction  or  a  multiple. 

The  failure  of  the  source  depths  to  beidentical  is  like’y  to  be 
the  most  common  and  significant  source  of  error  in  deconvolution  by  this 
method.  Figures  68-89  illustrate  the  effects  of  source  depth  differences 
on  the  deconvolution  of  a  one  bottom  bounce  signal  from  a  two  bottom  bounce 
signal.  Figure  68  shows  the  input  impulse  response  to  be  used  in  these 
synthetics.  Figure  69  shows  the  deconvolution  using  identical  91.4  m 
source  waveforms.  The  source  waveforms  and  impulse  responses  were  cal¬ 
culated  by  the  methods  discussed  in  Section  1.5.  Figure  70  shows  that  the 
impulse  response  can  be  recovered  arbitrarily  well  for  identical  source 
waveforms.  This  is  true  regardless  of  the  source  depth  used. 

Figures  71-80  show  the  effect  of  deconvolving  a  1  BB  from  a  2  BB 
waveform  with  244.1  m  mean  source  depths  in  which  the  1  BB  source  waveform 
represents  1,  3,  5,  7,  and  9%  shallower  depths  than  the  2  BB  source 
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FIGURE  68 

INPUT  IMPULSE  RESPONSE  FOR  SYNTHETIC  DECONVOLUTION 


J  \ 

INPUT  IMPULSE  RESPONSE 

91.4  m  SOURCE  WAVEFORM 


1  8B  SYNTHESIZED  WAVEFORM 


2  BB  SYNTHESIZED  WAVEFORM 


- A - - - - 

RECOVERED  IMPULSIVE  RESPONSE 

f - i - 1 - 1 - 1 - 1 - 1 

0  100  200  300 

TIME  —  msec 


FIGURE  69 
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FIGURE  72 
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FIGURE  78 
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FIGURE  79 
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FIGURE  80 
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FIGURE  82 

+7%  SOURCE  DEPTH  ERROR 

(244.1  m  SOURCE;  ALL  PASS) 
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FIGURE  83 

+7%  SOURCE  DEPTH  ERROR 

(244.1  m  SOURCE;  0-350  Hz) 
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waveform.  It  is  seen  for  all-pass  deconvolution  of  sources  with  mean  depths 
at  244  m  that  the  deconvolution  cannot  tolerate  greater  than  about  3%  source 
depth  differences. 


However,  as  shown  in  Figs.  81-85,  bandpass  filtering  reduces  the 
noise  level  at  the  cost  of  resolution.  Figures  81  and  82  show  an  all-pass 
deconvolution  of  waveforms  with  mean  source  depths  of  244.1  m  in  which  the 
1  BB  source  waveform  represents  7%  deeper  detonation  than  the  2  BB  source 
waveform.  It  is  seen  that  without  a  priori  knowledge  of  the  impuse  response 
the  specular  reflection  is  not  detectable.  Figures  83-85  show  bandpasses 
of  0-350  Hz,  0-120  Hz,  and  10-50  Hz,  respectively,  and  it  is  seen  that  the 
reflected  event  is  progressively  easier  to  detect  although  the  resolution 
decreases. 

Figures  86-89  show  the  effects  of  source  depth  differences  of  3t 
for  mean  depths  of  91.4  m  and  18.3  ®.  It  is  seen  that  3%  source  depth 
differences  make  the  specular  reflection  almost  undetectable. 

Finally,  Fig.  90  shows  an  example  of  the  deconvolution  of  five 
bottom  bounce  waveforms  from  six  bottom  bounce  waveforms.  The  same  multiple 
bottom  bounce  data  as  shown  in  previous  sections  of  this  report  were  used. 
The  experiment  geometry  was  a  244.1  m  source  and  a  256  m  receiver.  A  flat 
bandpass  was  applied  from  40-160  Hz.  This  clearly  identified  the  refracted 
arrival  and  a  reflected  arrival,  although  it  is  ambiguous  whether  this  is 
the  specular  reflection  or  a  subbottom  reflection.  The  delay  between  these 
two  components  is  the  reflection-refraction  delay  recovered  by  the  correla¬ 
tion  techniques.  The  relative  advance  of  the  specular  reflection  expected  a 
priori  for  this  turbidite  province  is  also  shown  in  the  figure.  The  data 
not  clearly  support  identification  of  an  arrival  at  these  times.  Energy 
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trailing  the  refracted  arrival  is  tentatively  identified  as  due  to  a  deeper 
refraction  and  multiple  reflection  events. 
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V 


CHAPTER  5 


GEOLOGIC  INFERENCES 

In  the  past  sections  the  correlation  techniques  and  the  impulse 
response  estimation  technique  imply  that  two  components  dominate  the  bot¬ 
tom  return,  a  refracted  component  and  a  reflected  component.  The  reflected 
component  may  be  either  the  specular  reflection  or  a  subbottom  reflection. 
If  it  is  a  subbottom  reflection  then  the  geoacoustic  structure  In  this 
province  is  significantly  different  than  that  predicted  a  priori  on  the 
basis  of  other  seismic  evidence  and  coring  data  (Hamilton  and  Bachman, 
personal  communication).  On  the  other  hand,  if  the  reflected  component  is 
the  specular  reflection,  it  will  be  shown  that  the  velocity  profile  for  the 
upper  few  hundred  meters  in  the  subbottom  is  significantly  greater  than 
that  predicted  a  priori.  Of  course  the  correlation  techniques  and  the 
impulse  response  estimates  indicated  that  it  is  likely  that  there  are  other 
components  of  the  ocean  bottom  transfer  function  besides  the  single 
reflection  and  the  refracted  component.  Although  these  other  arrivals  are 
ignored  in  obtaining  estimates  of  velocity  and  attenuation  as  a  function 
of  depth,  it  will  be  seen  that  the  results  imply  that  there  are  more  than 
two  components. 

The  technique  used  here  to  determine  the  velocity  profile  from 
the  time  delay  between  the  reflected  and  refracted  components  is  due  to 
Focke  (1981)  and  makes  use  of  a  computer  program  written  by  Smith  (1981). 
The  method  is  essentially  one  of  ray  tracing  the  reflected  and  refracted 
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paths  and  recognizing  that  the  ray  parameter  is  the  reciprocal  of  the 
velocity  at  the  turning  point  of  the  refracted  arrival,  i.e., 


(5.1) 


where  Cy  is  the  vertex  velocity,  r  is  the  range,  and  t  is  the  time. 

In  the  present  study,  since  the  arrival  time  of  the  refracted 
path  is  not  known,  but  rather  the  time  difference  between  the  reflected  and 
refracted  components,  Eq.  (5.1)  becomes 
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where  <t>  is  the  reflection  angle,  $w  is  the  refraction  angle,  and  Cw  is  the 
water  sound  speed  at  the  bottom.  The  travel  time  differences  were  deter¬ 
mined  from  the  correlation  techniques  of  Section  3.5  using,  in  particular, 
Fig.  56. 

The  refraction  angle  is  then  found  using 


oos(0w)  =  Cw  +  cos  (d>) 

The  range  through  the  sediment  for  the  refracting  ray  will  be  the  differ¬ 
ence  between  the  ranges  for  reflected  and  refracted  paths  through  the  water 
determined  by  ray  tracing. 


r(0w)  =  R (<P)  -  R(0m) 


152 


The  travel  time  through  the  sediment  is  determined  in  the  same  manner. 

r(6  )  =  At(<*>) +  T(<») -TvflJ  . 
w  w 

Since  the  time  delay  in  the  water  column  can  be  calculated  for  any  given 
angle,  knowledge  of  the  time  difference  between  the  reflected  and 
refracted  arrivals  allow  for  the  sound  speed  to  be  determined  for  a  range  r 
associated  with  an  angle  $ related  to  a  depth  h, 

r  sintf 

w 

k  2(cos0  +1) 
w 

The  sound  speed  at  this  depth  can  be  calculated  from  Eq.  (5.1). 

The  result  of  this  procedure  for  the  travel  time  differences 
between  the  reflected  ana  refracted  components  indicated  in  previous  sec¬ 
tions  is  shown  in  Fig.  91.  It  is  seen  that  the  sound  speeds  indicated  from 
the  multiple  bottom  bounce  signals  is  greater  than  that  predicted  a  priori 
for  the  upper  200  m  of  sediment.  This  sound  speed  profile  was  constrained 
at  the  water-sediment  interface  to  be  the  sound  speed  determined  a  priori 
from  coring  data  in  this  region.  Without  this  constraint  the  data  indi¬ 
cated  a  sound  speed  at  the  water  sediment  interface  of  about  1560  m/sec, 
which  is  a  sound  speed  ratio  at  the  water-sediment  interface  of  1.03  rather 
than  1.00.  Regardless  of  which  surface  value  is  used  the  sound  speed 
profiles  determined  from  the  time  delays  of  the  multiple  bottom  bounce  data 
merge  after  about  40  m  into  the  subbottom. 

The  velocity  profile  indicated  in  Fig.  91  for  the  multiple 

bottom  bounce  data  can  be  interpreted  in  terms  of  three  layers  with  linear 
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gradients  that  vary  continuously  across  the  interfaces.  The  first  layer 
varies  from  1515  m/sec  at  the  surface  to  1640  m/sec  at  60  m,  a  gradient  of 
2.08  sec-1.  The  third  layer  varies  from  1715  m/sec  at  110  m  to  1800  m/sec 
at  185  m,  a  gradient  of  1.13  sec-1. 

This  certainly  is  not  the  only  interpretation  of  the  variation 
of  sound  speed  as  a  function  of  depth  from  this  profile,  particularly  in 
the  upper  60  m,  but  it  is  a  fairly  simple  model  and  matches  most  of  the 
curve.  Using  this  velocity  profile,  an  attenuation  profile  for  the  turbi- 
dite  province  was  calculated  from  the  bottom  loss  data  using  the  method  of 
Mitchell  and  Focke  (1980). 

Figure  92  shows  the  attenuation  profile  calculated.  It  is  seen 
that  the  attenuation  varies  from  0.050  at  the  surface  to  about  0.015  at  a 
depth  of  about  160  m.  The  darker  line  represents  the  interpreted  attenua¬ 
tion  curve  for  the  attenuation  values  calculated  for  each  of  five  frequency 
bins  of  25,  50,  100,  200,  and  400  Hz.  At  shallow  depths  the  interpreted 
curve  should  be  biased  with  respect  to  the  higher  frequency  bins  since  it 
is  difficult  to  measure  the  low  loss  involved  at  the  lower  frequencies.  At 
the  greater  depths  the  interpreted  curve  should  be  biased  with  respect  to 
the  lower  frequency  bins  since  the  higher  frequencies  will  be  attenuated  to 
the  noise  level.  For  all  frequency  bins  other  than  400  Hz  the  attenuation 
curves  are  close  together  below  50  m.  Above  50  m  there  are  large  varia¬ 
tions  of  attenuation  as  a  function  of  depth.  This  is  interpreted  as 
evidence  of  unmodeled  reflectivity  in  the  upper  50  m  of  sediment.  Periodic 
occurrence  of  reflectors  in  the  subbottom  would  act  as  a  high  pass  filter 
of  the  ocean  bottom  transfer  function  while  reflections  due  to  sharp,  but 
continuous,  changes  in  impedance  with  depth  act  as  low  pass  filters. 
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Table  IV  shows  a  geoacoustic  model  derived  from  the  -elocity  and 
attenuation  profiles  of  Figs.  91  and  92.  The  densities  were  calculated 
from  the  velocity-density  relationships  of  Hamilton  for  terrigenous  sedi¬ 
ments.  Below  180  m  the  geoacoustic  model  was  found  by  connecting  a  line 
between  the  geoacoustic  parameters  at  this  depth  with  the  parameters  sug¬ 
gested  a  priori  at  1231  m  for  this  province.  Presumably,  there  is  a 
reflector  at  this  depth  identified  in  other  seismic  work. 

Table  V  shows  the  geoacoustic  model  derived  a  priori  for  this 
province.  It  is  seen  that  this  model  suggests  a  gradient  at  the  surface  of 
about  1.15  sec”1  and  attenuation  values  that  are  almost  an  order  of  magni¬ 
tude  greater  than  those  predicted  from  the  multiple  bottom  bounce  data. 
Note  also  that  no  subbottom  reflectors  are  identified  above  1231  m. 

In  the  rest  of  this  section  autocorrelograms  and  power  cepstra 
will  be  calculated  synthetically  from  the  geoacoustic  models  shown  in 
Tables  IV  and  V.  Comparing  the  synthetic  results  to  the  data  autocorrelo¬ 
grams  and  power  cepstra  lead  to  the  expected  conclusion  that  the  synthetics 
calculated  from  the  geoacoustic  model  derived  from  the  multiple  bottom 
bounce  data  agree  with  data  correlations  better  than  those  derived  using 
the  a  priori  model.  However,  it  is  also  found  that  there  are  significant 
differences  between  data  correlations  and  the  synthetic  correlations 
calculated  with  the  new  geoacoustic  model. 

Figure  93  shows  a  synthetic  autocorrelogram  calculated  from  the 
a  priori  geoacoustic  model.  The  reflection-refraction  component  and  the 
ghost  components  of  the  acoustic  signals  can  be  observed  on  the 
autocorrelogram  as  in  the  data  autocorrelogram  shown  in  Fig.  55.  However, 
there  are  some  significant  differences  between  the  autocorrelograms. 
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TABLE  IV 


NEW  GEOACOUSTIC  MODEL  DERIVED 
PROM  MULTIPLE  BOTTOM  BOUNCE  ACOUSTIC  DATA 


Depth 

(a) 

Velocity 

(n/sec) 

Attenuation 

(dB/m-kHz) 

Density 

(R/cm3) 

o- 

1514.5 

1 .04306 

0+ 

1515 

0.048 

1.59 

60 

1640 

0.020 

1.67 

110 

1715 

0.017 

1.76 

185 

1800 

0.015 

1.85 

1231” 

2443 

0.02 

2.26 

1231* 

2565 

0.02 

2.30 

172^* 

2900 

0.02 

2.40 

2217 

3735 

0.02 

2.49 

2217 

3400 

0.02 

2.52 

2270 

3685 

0.02 

2.57 

3820- 

3975 

0.02 

2.61 

3820+ 

4600 

0.02 

2.50 
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TABLE  V 

GEOACOUSTIC  MODEL  SUGGESTED  A  PRIORI 

(Shear  Waves  Ignored) 

Velocity 

Attenuation 

Density 

(m/see) 

(dB/m-kHz) 

(g/c»3) 

10 

1.59 

11 

1.67 

12 

1.80 

13 

1.91 

i4 

2.02 

14 

2.10 

12 

2.13 

1 

2.15 

0 

2.18 

18 

2.20 

>7 

2.22 
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2.24 

15 

2.25 

15 

2.26 

15 
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FIGURE  93 

SYNTHETIC  AUTOCORRELOQRAM  USING  A  PRIORI  GEOACOUSTIC 
MODEL  SHOWN  IN  TABLE  V 
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First,  the  synthetic  autocorrelogram  is  much  more  dominated  by  lower 
frequencies  than  the  data  autocorrelogram.  This  is  clear  evidence  that  the 
attenuations  used  in  the  a  priori  model  are  too  high.  In  addition, 
although  the  delays  of  the  ghosted  components  are  nearly  the  same  for  the 
two  autocorrelograms ,  the  reflection-refraction  delay  is  considerably 
different.  This  indicates  that  the  plane  wave  assumption  made  for  the 
calculation  of  the  ghost  delays  in  the  synthetic  data  is  approximately 
valid,  but  that  the  velocity  profile  in  the  sediment  is  inaccurate.  Thus, 
the  expected  conclusion  is  reached  that  the  velocity  and  attenuation  pro¬ 
files  calculated  for  this  region  a  priori  are  inconsistent  with  the  data. 

Figures  94  and  95  show  synthetic  power  cepstra  using  the  a 
priori  model.  It  is  seen  that  the  ghosts  and  reflection-refraction  delays 
can  be  detected  and  that  whitening  the  power  spectra  before  calculating  the 
power  cepstra  has  little  effect  on  the  synthetic  power  cepstrum  plots.  The 
amplitude  variation  of  those  components  as  a  function  of  grazing  angle  is 
dependent  somewhat  on  attenuation  along  the  refracted  path  as  evident  in 
the  autocorrelogram,  Fig.  93«  On  the  other  hand,  some  of  these  variations 
such  as  those  that  go  from  high  amplitude  to  low  amplitude  back  to  high 
amplitude  again  are  due  to  the  coupling  between  quefrency  components.  This 
illustrates  the  major  disadvantage  of  the  power  cepstrum,  which  is  espec¬ 
ially  significant  as  the  number  of  components  in  the  ocean  bottom  transfer 
function  becomes  large  (beating  patterns  become  very  complicated  and 
amplitude  variations  become  hard  to  predict). 

On  the  other  hand,  the  major  advantage  of  the  power  cepstrum  over 
the  autocorrelation  is  that  a  single  spike  represents  the  time  delay 
between  components  with  an  amplitude  determined  by  the  average  dynamic 
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range  of  the  ripple.  For  some  signals  this  makes  the  time  delay  easier  to 
detect  using  the  power  cepstrum  than  using  the  autocorrelation,  which 
gives  a  pulse  at  the  time  delay  between  components  with  a  power  spectrum 
determined  by  the  crosspower  spectrum  of  the  impulse  response  components. 
Thus,  if  the  correlated  energy  of  the  impulse  response  is  dominated  by  low 
frequencies  the  delay  will  be  hard  to  resolve  in  the  autocorrelation.  This 
effect  can  be  seen  by  comparing  Figs.  93  and  94  where  the  reflection- 
refraction  delay  is  detected  easier  on  the  power  cepstrum  than  the  auto¬ 
correlation  for  angles  less  than  16°. 

Comparing  the  synthetic  power  cepstra  calculated  using  the  a 
priori  model  (Figs.  94  and  95),  with  the  data  power  cepstra  (Figs.  60-63), 
it  is  seen  that  there  are  some  significant  differences.  First,  there  are 
many  more  cepstral  peaks  in  the  data  than  in  the  synthetics.  This  is 
believed  to  be  due  primarily  to  shallow  reflections,  for  instance  in  the 
upper  50  m,  although  it  may  also  be  due  to  the  effects  of  beating  between 
quefrencies  and  additive  random  noise,  which  may  be  colored. 

Second,  many  of  the  peaks  found  both  in  the  data  and  in  the 
synthetics  are  of  notably  different  amplitudes  and  shapes.  For  instance, 
the  first  bubble  pulse  peak  in  the  data  is  very  large,  whereas  the  first 
bubble  pulse  peak  is  barely  visible  in  the  synthetics.  Also  the  large  peak 
at  very  low  quefrencies  is  much  broader  in  the  synthetics  than  in  the  data. 
Finally,  the  peak  associated  with  the  reflection-refraction  delay  is  not 
bimodal  nearly  as  soon  in  the  data  as  in  the  synthetics.  All  of  these 
effects  are  believed  to  be  due  primarily  to  the  fact  that  the  attenuation 
values  in  the  a  priori  model  are  too  high.  The  fact  that  the  ghost  delays 
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are  more  poorly  defined  in  the  data  is  believed  to  be  due  to  greater 
subbottom  reflectivity  and  an  attendant  increase  in  quefrency  beating. 

Figures  96-98  show  the  synthetic  correlations  using  the  new 
geoacoustic  model  shown  in  Table  IV.  There  are  some  similarities  and  some 
important  differences. 

Comparing  the  autocorrelograms,  Figs.  55  and  96,  it  is  seen  that 
the  main  components  of  energy  in  the  autocorrelogram  of  the  data  are  repro¬ 
duced  in  the  synthetic  autocorrelogram.  The  quiet  regions  on  the  auto¬ 
correlogram  of  the  data  are  quiet  on  the  synthetic  autocorrelogram.  This 
shows  that  the  source-receiver  geometry  is  well  modeled  in  the  data  and  the 
main  components  of  the  ocean  bottom  transfer  function  are  modeled. 

There  are,  however,  some  significant  differences  between  the 
autocorrelograms.  First  the  autocorrelations  of  the  data  at  angles  less 
than  about  26°  have  significantly  greater  high  frequency  content  than  the 
synthetic  autocorrelations  at  the  same  angles.  This  indicates  that  the 
attenuation  values  for  the  new  geoacoustic  model  are  too  high  at  the 
shallow  depths. 

The  occurrence  of  high  frequencies  at  low  angles  for  the  data 
autocorrelogram  may  be  due  to  reflections  at  shallow  depths.  The  occur¬ 
rence  of  reflections  has  the  effect  of  a  high  pass  filter  since  low  fre¬ 
quencies  do  not  see  reflectors  spaced  considerably  closer  than  their 
wavelengths.  Reflectors,  if  they  ocurr,  must  be  modeled  if  an  accurate 
attenuation  profile  is  to  be  determined. 

Another  difference  between  the  data  and  synthetic 
autocorrelograms  is  the  time  delay  between  the  reflected  and  refracted 
components  at  angles  greater  than  about  30°.  This  is  due  to  the  fact  that 
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FIGURE  96 

SYNTHETIC  AUTOCORRELOGRAM  USING  NEW  GEOACOUSTIC  MODEL 
PRESENTED  IN  TABLE  IV  DERIVED  FROM  DATA  CORRELOGRAMS 
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FIGURE  97 

SYNTHETIC  POWER  CEPSTRA  USING  NEW  GEOACOUSTIC  MODEL 
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FIGURE  98 

8YNTHETIC  POWER  CEP8TRUMS  USING  NEW  GEOACOUSTIC  MODEL 

WITH  WHITENED  POWER  SPECTRA  ARL  UT 
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the  reflection-refraction  delay  for  the  synthetics  was  calculated  on  the 
basis  of  a  plane  wave  assumption  which  breaks  down  at  the  higher  angles. 
In  addition,  for  real  data  with  spherical  wavefronts,  there  is  a  minimum 
range  or  maximum  grazing  angle  for  which  a  refracted  arrival  will  occur 
(Spofford,  1980).  Thus,  the  reflection-refraction  delay  of  the  synthetic 
autocorrelogram  can  be  detected  for  angles  up  to  almost  40°,  while  the 
reflection-refraction  component  of  the  data  autocorrelogram  cannot  be  seen 
beyond  about  32°.  The  failure  of  the  reflection-refraction  component  to 
occur  could  be  due  to  the  fact  that  no  refraction  component  occurs  beyond 
that  angle  or  to  the  fact  that  it  is  attenuated  below  the  noise  level. 

Finally,  in  this  section,  the  synthetic  power  cepstra  using  the 
new  geoacoustic  model  (Figs.  97  and  98)  will  be  compared  to  the  data  power 
cepstra.  The  amplitudes  of  the  peaks  of  the  main  cepstral  components  of 
the  synthetic  power  cepstra  are  in  much  better  agreement  with  the  data 
than  the  synthetic  power  cepstra  derived  using  the  a  priori  model.  For 
instance,  the  amplitude  of  the  bubble  pulses,  the  relative  amplitude  of  the 
ghost  components,  and  the  amplitude  and  shape  of  the  reflection-refraction 
component  are  more  nearly  the  same  than  they  are  for  the  a  priori  model. 
Even  the  relative  highs  and  lows  in  the  amplitudes  of  the  ghost  components 
(due  to  beating)  is  duplicated  in  the  new  model,  especially  at  angles 
between  about  10  and  18°.  This  agreement  is  due  to  the  lower  attenuations 
assumed  in  the  new  model. 

The  differences  between  the  synthetic  power  cepstra  and  the  data 
are  primarily  that  there  are  many  cepstral  peaks  in  the  data  at  delays 
where  there  are  none  in  the  synthetics.  This  is  believed  to  be  due  to 


reflections  and  their  multiples  and  to  beating  between  quefrencies 
associated  with  these  components. 

In  summary,  the  new  geoacoustic  model  is  an  improvement  over  the 
previous  model  for  this  region  since  it  improves  the  match  of  the  delay 
between  the  two  main  components  of  the  ocean  bottom  impulse  response  and 
the  frequency  content  of  the  autocorrelations  as  a  function  of  grazing 
angle.  However,  it  is  still  deficient  in  two  areas.  First,  the  attenua¬ 
tion  values  may  be  still  lower  at  the  shallow  depths  than  those  calculated 
earlier  in  this  section  and,  second,  there  may  be  one  or  more  subbottom 
reflectors  at  the  shallow  depths.  Accurate  determination  of  the  attenua¬ 
tion  coefficient  is  dependent  on  knowledge  of  the  occurrence  of 
reflectors. 
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CHAPTER  6 


CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  Conclusions 

In  the  first  chapter  of  this  thesis  an  overview  of  bottom  loss 
collection  was  given  and  the  state-of-the-art  knowledge  of  geoacoustic 
parameters  for  different  ocean  environments  was  presented.  It  was  empha¬ 
sized  that  even  in  the  deep  ocean  the  physical  and  geologic  processes  are 
complicated  and  variable,  thus  giving  considerable  variability  to 
geoacoustic  structure.  This  variability  is  seen  by  examining  physical 
properties  of  geologic  cores  taken  in  relatively  close  proximity,  by  com¬ 
paring  bottom  loss  curves  in  seemingly  similar  ocean  environments,  and  by 
comparing  deconvolved  acoustic  measurements,  which  reveal  variable  and 
quite  significant  subbottom  reflectivity  structure  (Hastrup,  1970; 
Herstein  et  al,  1979;  Santaniello  et  al,  1979;  Chapman,  1980;  Tyce  et  al, 
1980;  Dicus,  1981).  Topographic  effects  such  as  the  depth  and  nature  of 
the  basement  also  greatly  affect  bottom  loss.  In  shallow  water  it  would  be 
expected  that  topographic  effects  would  also  greatly  influence  bottom 
loss.  It  would  be  expected  that  the  geoacoustic  structure  would  be  even 
more  complex  and  variable  in  shallow  water  than  in  deep  water. 

In  Chapter  2  it  was  shown  that  multiple  bottom  bounce  acoustic 
signals  can  be  described  by  a  series  of  linear  filters  which  include  bubble 
pulses  and  ghosting  effects.  The  effects  of  ghosts  and  bubble  pulses  can 
often  be  identified  on  the  acoustic  data  and,  in  general,  they  severely 


distort  the  ocean  bottom  impulse  response.  For  this  reason  it  is  difficult 
to  extract  information  on  the  geoacoustic  properties  of  sediments  in  the 
ocean  bottom  from  the  unprocessed  acoustic  signals.  It  has  been  shown  that 
the  ghosts  and  bubble  pulses  can  be  adequately  identified  and  that  the 
nominal  source  and  receiver  depths  can  be  checked.  It  was  found  that,  in 
the  data  discussed  in  this  report,  the  actual  receiver  depth  (~255  m)  was 
considerably  less  than  the  nominal  value  (498  m). 

The  two  correlation  techniques  investigated  in  this  work,  the 
autocorrelation  and  power  eepstrum,  greatly  facilitate  analysis  of  the 
acoustic  signals.  The  ghost  delays  are  detectable  in  many  cases  even  with 
experimental  geometries  that  were  not  favorable  for  their  detection,  and 
an  estimate  of  the  variance  of  the  source  depth  is  obtainable  from  the 
bubble  pulse  delay.  In  addition,  inferences  into  the  geoacoustic  struc¬ 
ture  are  obtainable  from  the  correlations;  a  delay  associated  with  a 
reflected  and  the  refracted  component  is  detectable.  Initially  the  auto¬ 
correlation  is  preferable  to  the  power  eepstrum  since  it  has  the 
superposition  property  with  respect  to  convolution,  and  thus  is  easier  to 
interpret  than  the  power  eepstrum.  Coherent  noise  such  as  ghosts  and 
bubble  pulses  are  more  easily  recognized,  the  amplitudes  of  events  are  less 
variable,  and  the  shapes  of  events  are  more  interpretable  in  terms  of  the 
phase  of  the  time  arrivals.  All  of  these  characteristics  are  due  to  the 
fact  that  the  power  eepstrum  has  a  nonlinear,  logarithmic  operation  and 
therefore  there  is  coupling  between  cepstral  components.  This  is  what 
causes  the  unpredictable  variability  of  the  amplitude  of  the  cepstral 
components.  Also,  the  shapes  of  events  are  not  duplicated  on  the  power 
oepstrum  causing  some  difficulty  in  recognizing  signal 
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components,  and  the  power  cepstrum  is  significantly  degraded  by 
uncorrelated  noise. 

It  was  found,  however,  that  for  much  synthetic  and  real  data  the 
power  cepstrum  was  as  good  and  better  than  the  autocorrelation  for 
measuring  time  delays  between  components  of  the  ocean  bottom  impulse 
response.  This  is  due  to  the  fact  that  the  power  cepstrum  gives  a  single 
spike  at  the  time  delay  between  components  with  an  amplitude  determined  by 
the  dynamic  range  of  the  ripple,  whereas  the  amplitude  of  the  autocorrela¬ 
tion  is  determined  by  the  crosspower  spectrum  of  the  components,  which  may 
be  highly  dominated  by  low  frequencies  leading  to  poor  resolution.  Thus, 
especially  at  low  grazing  angles  where  the  S/N  ratio  is  high  and  the  time 
delay  between  components  is  small,  the  power  cepstrum  may  be  preferable. 

The  method  of  deconvolving  an  (n-1)  bottom  bounce  signal  from  an 
(n)  bottom  bounce  signal  is  appealing  due  to  the  ease  of  implementation. 
It  has  been  shown  to  be  effective  in  estimating  the  main  components  of  the 
impulse  response  over  a  fairly  limited  bandwidth.  Criteria  that  are 
crucial  for  this  method  of  deconvolution  are  nearly  identical  source  and 
receiver  depths,  and  sufficient  number  of  bottom  bounces  that  the  ghost 
angles  are  nearly  identical  between  the  (n)  and  (n-1)  bottom  bounce  wave¬ 
forms. 

In  the  turbidite  region  examined  in  this  study  there  appear  to  be 
two  main  components  of  the  impulse  response— a  reflected  event  and  a 
refracted  event.  It  is  not  known  whether  the  reflector  is  the  specular 
reflector  or  a  subbottom  reflector,  but  in  either  case  the  results  imply 
that  the  geoacoustic  structure  is  significantly  different  than  that  pre¬ 
viously  predicted.  If  the  reflector  is  the  specular  reflector  then  the 
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sound  speed  increases  with  depth  much  more  rapidly  than  previously 
predicted;  otherwise,  there  is  a  significant  discontinuity  or  a  sharp 
increase  of  impedance  with  depth  that  was  previously  not  modeled. 

In  addition  there  is  evidence  from  the  power  cepstrum  and  18  m 
source  autocorrelations  that  there  is  significant  subbottom  reflectivity. 
The  attenuation  profiles  at  different  frequencies  seem  to  confirm  signifi¬ 
cant  reflectivity  in  the  upper  40  m  of  the  subbottom.  This  is  consistent 
with  the  recent  work  of  Tyce  et  al.  ( 1 980)  that  indicates  strong  reflectors 
of  the  order  of  the  specular  reflector  up  to  40  m  below  the  water-sediment 
interface.  The  S/N  ratio  of  the  deconvolved  traces  is  not  sufficient  to 
identify  with  confidence  the  reflectivity  implied  by  the  autocorrelations 
and  power  cepstra. 

Assuming  that  the  identified  reflector  is  the  specular 
reflector,  the  gradient  at  the  surface  is  higher  (g~2.1)  than  that  gen¬ 
erally  assumed  (g~1.2),  although  gradients  at  the  surface  greater  than 
1.2  see”1  have  been  measured  recently,  e.g.,  1.6  sec-1  (Shirley,  1978; 
Tucholke,  1980)  and  1.8  sec”1  (Dorman,  1981).  For  the  data  shown  here  the 
gradient  decreases  to  about  1.13  sec”1  at  about  180  m  depth.  This  gives  an 
average  gradient  of  about  1.5  sec”1  for  the  upper  100  m.  This  is  not 
greatly  different  from  the  value  of  1.35  sec”1  recently  reported  by 
Dinapoli  nor  those  used  by  Spofford  to  fit  bottom  loss  curves  of  gQ  s  1.7 
see”1  and  g(500)  a  1.34  sec”1. 

(Bottom  Interacting  Ocean  Acoustics.  1980). 

The  attenuation  values  which  were  derived  to  model  the 
autocorrelations  of  the  data  Indicate  that  k  ranges  from  about  0.05  to 
0.01  dB/m-kHz,  which  is  within  the  range  measured  by  Mitchell  and  Fooke 
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(I960),  although  these  values  are  an  order  of  magnitude  less  than  those 
measured  by  some  others  (Hamilton,  1976,  1980);  Jacobson,  1981).  Chapman 
(1980)  showed  that  low  frequency  bottom  loss  measurements  in  the  Tufts 
abyssal  plain  implied  an  attenuation  coefficient  of  0.002  dB/m  at  160  Hz 
(k~0.013)  and  Spofford  (1980)  required  an  attenuation  coefficient  which 
varied  from  0.019  at  the  surface  to  0.119  at  500  m  to  fit  bottom  loss 
curves.  Thus,  it  is  concluded  that  the  attenuation  and  sound  speed  pro¬ 
files  derived  from  the  bottom  loss  data  are  not  unreasonable  with  respect 
to  the  state-of-the-art  knowledge  of  these  parameters. 

It  is  concluded  that  the  autocorrelation,  power  cepstrum,  and 
multiple  bottom  bounce  deconvolution  techniques  are  valuable  tools  for  the 
extraction  of  geoacoustic  information  from  multiple  bottom  bounce  signals. 
6.2  Recommendations 

There  are  several  recommendations  that  can  be  made  on  the  basis 
of  this  work.  First,  it  is  recommended  that  the  bottom  loss  curves  be 
recalculated  for  this  site  using  the  empirically  derived  receiver  depths. 
The  difference  between  the  nominal  and  empirical  ghost  time  delays  will  be 
particularly  significant  at  the  higher  grazing  angles  and  this  error  will 
be  reflected  in  the  scalloping  in  the  frequency  domain. 

Second,  it  is  recommended  that  study  of  the  correlation 
techniques  continue.  A  particularly  interesting  study  would  be  to  compare 
autocorrelograms  over  a  range  of  passbands  for  different  ocean  basins 
where  multiple  bottom  bounce  autocorrelations  of  the  same  source  and 
receiver  depths  and  bottom  grazing  angles  are  summed.  Muting  along  the 
source  or  receiver  ghost  delay  before  autocorrelating  would  be  worthwhile 
in  improving  the  resolution  of  correlated  events.  Study  of  power 
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cepstra  for  different  ocean  basins  would  also  be  of  value,  particularly  in 
examining  subbottom  reflectivity. 


Third,  improvement  of  the  deconvolution  technique  is  needed.  It 
would  be  worthwhile  to  apply  the  technique  of  Dicus  (1981),  especially  if  a 
good  estimate  of  the  source  spectrum  can  be  obtained  by  measuring  the 
source  waveform  without  interference  of  ghosts  or  bottom  interacting  sig¬ 
nals.  The  shot  depth  variation  can  be  monitored  by  the  correlation  tech¬ 
niques.  If  a  good  estimate  of  the  source  spectrum  is  not  available,  then 
significant  Improvement  of  the  impulse  response  estimate  might  be  obtained 
by  posing  the  problem  as  predictive  deconvolution  with  operator  length 
less  than  one. 

In  addition  a  few  recommendations  can  be  made  with  regard  to  the 
experimental  geometry.,  First,  deeper  sources  (>400  m)  would  greatly  faci¬ 
litate  study  of  the  acoustic  signals.  This  would  allow  for  accurate  muting 
of  the  data  beyond  the  source  ghost,  especially  if  time  separation  between 
Impulse  response  components  with  ghost  components  occurs.  Also,  the 
resolution  of  the  impulse  response  estimate  would  be  Increased  since  the 
time  delay  between  bubble  pulses  is  almost  inversely  proportional  to 
depth,  and  so  the  errors  of  the  bubble  pulse  delays  are  almost  inversely 
proportional  to  source  depth  errors.  Thus,  the  model  breakdown  noise  of 
the  deconvolution  procedure  will  be  restricted  to  higher  and  higher  fre¬ 
quencies  as  the  source  detonation  depth  increases. 

Alternatively,  the  deconvolution  procedure  could  be  improved 
greatly  by  substantial  reduction  of  the  variance  of  the  source  detonation 
depths.  This  could  be  done  by  improving  the  quality  of  the  pressure  trans¬ 
ducers.  A  cost  effectiveness  study  would  need  to  be  made  to  choose  between 


increasing  the  pressure  transducer  quality  or  increasing  the  source 
detonation  depth. 

If  the  desired  data  product  is  an  impulse  response  estimate 
versus  range  and  angle  then  the  source  should  be  detonated  in  even  incre¬ 
ments  of  range  assuming  that  the  water  depth  and  sound  speed  is  invariant 
with  range.  The  size  of  the  increment  is  dependent  on  the  water  depth  and 
sound  speed  and  the  desired  range  and  angle  sampling.  This  would  allow  for 
an  optimum  estimate  of  the  impulse  response  for  multiple  bottom  bounce 
signals  in  which  the  source  spectrum  is  not  known.  Of  course,  conventional 
techniques  could  be  used  if  the  source  spectrum  is  known  or  can  be  esti¬ 
mated  accurately  (Wood  et  al.,  1979;  Dicus,  1981).  The  deconvolved  trace 
would  be  an  estimate  of  the  multiple  bottom  bounce  impulse  response  with 
its  source  and  receiver  ghosts.  The  ghosts  could  be  removed  by  predictive 
deconvolution  if  they  are  not  time  separated  from  the  impulse  response 
components.  An  estimate  of  the  one  bottom  bounce  impulse  response  could  be 
attempted  by  taking  the  nth  root  of  the  spectrum  of  the  (n)  bounce  impulse 
response  estimate. 

Finally,  the  whole  problem  of  determining  the  factors  which 
influence  angle  dependent  bottom  loss  would  be  greatly  facilitated  by  use 
of  array  deta,  for  instance,  in  a  multiship  experiment  such  as  ESP 
(expanding  spread)  or  fixed,  large  offset  shooting.  Tau-p  processing,  for 
instance,  could  then  be  used  to  model  the  arrival  structure  as  a  function 
of  grazing  angle,  which  determines  frequency  dependent  bottom  loss  as  a 
function  of  grazing  angle. 
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APPENDIX  I 


All  of  the  travel  time  curves  calculated  for  the  overlays  in 
Chapter  2  make  use  of  the  plane  wave  assumption.  For  travel  paths  in  the 
water  column  this  is  a  very  good  approximation  to  reality.  On  the  other 
hand,  paths  through  the  sediment,  especially  at  high  grazing  angles,  are 
more  accurately  modeled  assuming  spherical  wavefronts.  In  this  section,  a 
comparison  of  travel  time  curves  based  on  the  plane  wave  assumption  are 
compared  to  a  spherical  wave  model  in  an  ocean  model  with  an  isovelocity 
water  layer. 

Figure  A-1  shows  the  travel  time  difference  curves  for  the 
impulse  response  of  a  six  bottom  be..  -*.6  waveform  with  a  gradient  of 
1.5  sec-1,  a  sound  speed  ratio  of  one,  and  a  water  depth  of  3600  m.  The 
solid  curves  are  the  travel  time  curves  for  the  plane  wave  assumption  and 
the  dashed  curves  are  derived  assuming  spherical  waves  with  source  and 
receiver  at  the  top  of  the  water  column.  It  is  seen  that  the  travel  time 
curves  agree  very  well  for  signals  which  have  interacted  with  the  bottom 
only  a  few  times  at  angles  less  than  25°.  For  grazing  angles  greater  than 
25°  and  for  signals  which  have  interacted  with  the  bottom  several  times, 
the  travel  time  curves  have  significant  differences. 

Figure  A-2  shows  a  phenomenon  for  spherical  waves  that  does  not 

occur  for  plane  waves;  namely,  at  very  low  angles  spherical  wave  travel 

refracted  through  the  sediment  arrive  sooner  than  reflected  paths.  For  the 

same  parameters  as  used  in  Fig.  A-1,  except  c  /c  =0.98,  it  is  seen  that  the 

s  w 

reflected  and  refracted  paths  arrive  at  the  same  time  at  10°.  This  angle 
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will  be  dependent  on  the  source-receiver  geometry  and  the  sound  speed 
structure.  With  the  plane  wave  assumption  the  refracted  arrival  will  never 
arrive  before  the  reflected  arrival.  For  spherical  waves  the  refracted 
arrival  may  arrive  before  the  reflected  arrival  at  both  low  angles  and  at 
high  angles,  dependent  strongly  on  source-receiver  geometry  and  the  sedi¬ 
ment  sound  speed  structure. 

Finally,  in  connection  with  spherical  wave  fronts  it  should  be 
mentioned  that  there  is  a  minimum  range  for  a  given  source-receiver  geo¬ 
metry  and  bottom  sound  speed  structure  at  which  a  refracted  arrival  will 
occur  (Spofford,  1980). 


APPENDIX  II 


It  is  desired  to  show  that  the  autocorrelation  has  the 
superposition  property  with  respect  to  convolution  and  that  the  power 
cepstrum  is  not.  Consider  the  autocorrelations  of  the  functions  w(t)  and 
y(t). 


Pw(t)  s  C  w(u)  w(u-t)du 

—  OO 


Pyy(t)  =  f  y(u)  y(u+t)du 

*oo 


pww<tJ  •  Pyy(t) 


■£  £ 


w(u)  w(u-x)du 


][£>'■■ 


)  y(u’+t-x)du' 


dx 


Set  u’  s  -uf  in  second  inner  integral;  then 


P  (t)  •  P  (t) 
ww  yy 


w(u)  w(u-x)  y(-u)  y(t-u-x)du  dx 


•£['•"£ 

a  J  (-1)  w(u)  y(t-u-x)du||y*  w(u-x)  y(-u)duj 


Set  u’  a  u-x  in  second  inner  integral. 


Pww(t)  *  Pyy(t) 


s  (-1)  ^  w(u)  y  (t-u-x)dujj^j*  w(u’)  yO-u'-xJdu'^dx 
=  y  (-l)^w(t-x)  *  y(t-x)J£w(-x)  *  y(-x)Jdx 


Finally  set  x*  a  t-x, 


w(x')  *  y(x* 

Let  w(t)*y(t)=z(t). 

Pww(t)  •  PyyCt)  =  /°°z(x‘)  z(x'-t)  dx»  =  P2Z(t) 

Thus,  the  convolution  of  the  autocorrelations  of  two  functions  is  equal  to 
the  autocorrelation  of  the  convolution  of  two  functions,  and  the  auto¬ 
correlation  has  the  superposition  property  with  respect  to  convolution. 

Taking  the  Fourier  transform  of  the  convolution  of  the 
autocorrelations  yields  a  product  of  two  power  spectrums.  The  logarithmic 
operation,  however,  does  not  have  the  superposition  property  with  respect 
to  multiplication;  therefore,  the  power  cepstrum  does  not  have  the 
superposition  property  with  respect  to  convolution. 


pw“)  •  V‘>  ’ll 


)J  jw(x’-t)  •  y(x»-t)J 
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